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ABSTRACT 

Context. LOFAR offers the unique capability of observing pulsars across the 10 - 240 MHz frequency range with a fractional band¬ 
width of roughly 50%. This spectral range is well suited for studying the frequency evolution of pulse profile morphology caused 
by both intrinsic and extrinsic effects such as changing emission altitude in the pulsar magnetosphere or scatter broadening by the 
interstellar medium, respectively. 

Aims. The magnitude of most of these effects increases rapidly towards low frequencies. LOFAR can thus address a number of open 
questions about the nature of radio pulsar emission and its propagation through the interstellar medium. 

Methods. We present the average pulse profiles of 100 pulsars observed in the two LOFAR frequency bands: high band (120- 
167 MHz, 100 profiles) and low band (15-62 MHz, 26 profiles). We compare them with Westerbork Synthesis Radio Telescope 
(WSRT) and Lovell Telescope observations at higher frequencies (350 and 1400 MHz) to study the profile evolution. The profiles 
were aligned in absolute phase by folding with a new set of timing solutions from the Lovell Telescope, which we present along with 
precise dispersion measures obtained with LOFAR. 

Results. We find that the profile evolution with decreasing radio frequency does not follow a specific trend; depending on the geometry 
of the pulsar, new components can enter into or be hidden from view. Nonetheless, in general our observations confirm the widening 
of pulsar profiles at low frequencies, as expected from radius-to-frequency mapping or birefringence theories. We offer this catalogue 
of low-frequency pulsar profiles in a user friendly way via the EPN Database of Pulsar Profiles*. 

Key words, stars: neutron - (stars) pulsars: general 


1. Introduction 

: cumulative (i.e. average) pulse profiles of radio pulsars are 
sum of hundreds to thousands of individual pulses, and are, 
loosely speaking, a unique signature of each pulsar (iLorimeil 
l2008l) . They are normally stable in their morphology and are 
reproducible, although several types of var iation have been ob- 


served both for non-recvcled pulsars (see 

Helfand etal.ll 19751 

Weisberg et al.lll989llRathnasree & RankinI 

19951 and Lvne et al.l 

2010|) and for millisecond pulsars (|Liu et 

al.ll2012h. For most 


pulsars, this cumulative pulse profile morphology often varies 
(sometimes drastically, sometimes subtly) as a function of ob¬ 
serving frequency because of a number of intrinsic effects (e.g. 
emission location in the pulsar magnetosphere) and extrinsic ef¬ 
fects (i.e. due to propagation in the interstellar medium; ISM), 

* http://www.epta.eu.org/epndb/ 


see for instance fCordesI (Il978l) . As Fig. [T] shows, pulse profile 
evolution can become increasingly evident at the lowest observ¬ 
ing frequencies (< 200 MHz). 

Mapping profile evolution over a wide range of frequencies 
can aid i n modelling th e pulsar radio emission mechanism itself 
(see e.g. lRankinl[l993l and related pa pers of the series) a nd con¬ 
straining properties of the ISM (e.g. iHassall et al.ll2012l) . Since 
many of the processes that affect the pulse shape strongly depend 
on the observing frequency, observations at low frequencies pro¬ 
vide valuable insights on them. The LOw-Frequency ARray 
(LOFAR) is the first telescope capable of observing nearly the 
entire radio spectrum in the 10 - 240 MH z frequency range, 
the lo west 4.5 octaves of the ‘ra dio window’ (Ivan Haarlem et aH 
12013l) and (IStappers et alJl20lTl) . 

Low-frequency pulsar observations have previously been 
conducted by a number of telescopes, (e.g. Gauribidanur Ra- 
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Fig. 1. Example of pulsar profile evolution for PSR B0950+08, from 
1400 MHz down to 30 MHz. It becomes more rapid at low frequencies. 
The bars on the left represent the intra-channel smearing due to uncor¬ 
rected DM delay within a channel at each frequency. The profiles were 
aligned using a timing ephemeris (see text for details). 


dio Telescope (GEETEE); lAsgekar & Deshpande 199^ Larg e 
Phased Array Radio Teles cope. Pusc hino: Kimnin_et^ (H99®, 
iMalov & Malofee^ 120101 : Arecibo: iHankins & Rankinl 20101 : 
Ukrainia n T-shaped Radio tele scope, second modification 
(UTR-2): IZakharenko et al.ll2013h . and simultaneous efforts are 
being undertaken by other groups in par allel (e.g. Long Wave¬ 
length Array ( LWA): IS tovall et^ 2015L Murchison Widefield 
Array (MWA): lTremblav et al.l 2015h . Nevertheless. LOEAR of¬ 
fers several advantages over the previous studies. Eirstly, the 
large bandwidth that can be recorded at any given time (48 MHz 


in 16-bit mode and 96 MHz in 8-bit mode) allows for continu¬ 
ous wide-hand studies of the pulse profile evolution, compared 
to studies using a number of widely separated, narrow bands 
(e.g. the 5 x 32 an d 20 x 32 kHz bands used at 102 MHz in 
iKuzmin et al.lll998l) . Secondly, LOEAR’s ability to track sources 
is also an advantage, as many pulses can be collected in a single 
observing session instead of having to combine s everal short ob¬ 
serva tions in the case of transit instruments (e.g. lizvekova et aP 
Il989h . This eliminates systematic errors in the profile that are 
due to imprecise alignment of the data from several observing 
sessions. Thirdly, LOEAR can achieve greater sensitivity by co¬ 
herently adding the signals received by individual stations, giv¬ 
ing a collecting area equivalent to the sum of the collecting area 
of all stations (up to 57() 00m^ at HBAs and 75200 m^ at LBAs, 
see IStappers et al.ll201 Ih . Einally, LOEAR offers excellent fre¬ 
quency and time resolution. This is necessary for dedispersing 
the data to resolve narrow features in the profile. LOEAR is also 
capable of coherently dedispersing the data, although that mode 
was not employed here. 

As mentioned above, there are two types of effects that LO¬ 
EAR will allow us to study with great precision. One of these are 
extrinsic effects related to the ISM. Specifically, the ISM causes 
scattering and dispersion. Mean scatter-broadening (assuming a 
Kolmogorov distribution of the turbulence in the ISM) scales 
with observing frequency as The scattering causes delays 
in the arrival time of the emission at Earth, which can be mod¬ 
elled as having an exponentially decreasing probability of being 
scattered back into our line-of-sight: this means that the inten¬ 
sity of the pulse is spread in an exponentially decreasing tail. 
Dispersion scales as and is mostly cor rected for by chan¬ 
nelisi ng and dedispersing the data (see e.g. iLorimer & Kramen 
l2004l) . Nonetheless, for filterbank (channelised) data some resid¬ 
ual dispersive smearing persists within each channel: 

foAT = 8.3 • DM—jUS, (1) 


where DM is the dispersion measure in cm“^ pc, Av is the chan¬ 
nel width in MHz, and v the central observing frequency in GHz. 
This does not significantly affect the profiles that we present here 
(at least not at frequencies above ~ 50 MHz) because the pulsars 
studied have low DMs and the data are chanellised in narrow fre¬ 
quency channels (see Sect. |2]i. Second-order effects in the ISM 
may also be present, but have yet to be confirmed. Eor instance, 
previous claims of ‘super-di spersion’, meaning a deviation from 
the scaling law (see e.g. iKuz’min et aDl200^ an d references 
therein), were not observed bv iHassall et al.l (l2012h . with an up¬ 
per limit of < 50 ns at a reference frequency of 1400 MHz. 

The second type of effects under investigation are those in¬ 
trinsic to the pulsar. One of the most well-known intrinsic ef¬ 
fects are pulse broadening at lo w frequencies, which has b een 
observed in many pu lsars (e.g. IHankins & Ricketl Il986l and 
iMitra & Rankinl l2002h . while others show no evidence of this 
(e.g. Hassall et al. 2012). One of the theori es explaining this 
effect is radius-to-frequency mapping (REM. ICordeslflOT^ : it 
postulates that the origin of the radio emission in the pulsar’s 
magnetosphere increases in altitude above the magnetic poles 
towards lower frequencies. REM predicts that the pulse profile 
will increase in width towards lower observing frequency, since 
emission will be directed tangentially to the diverging magnetic 
field lines of the magnetosph ere that corotates with the pulsar. 
An alternative interpretation (iMcKinnonI 1 1 997h proposes bire¬ 
fringence of the plasma above the polar caps as the cause for 
broadening: the two propagation modes split at low frequencies, 
causing the broadening, while they stay closer together at high 
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frequencies, causing depolarisation (thi s is investigated in the 
LOFAR work on pulsar polarisation, see lNoutsos et alfcOlSh . 

In this paper, we present the average pulse profiles of 100 
pulsars observed in two LOFAR frequency ranges: high band 
(119-167 MHz, 100 profiles) and low band (15-62 MHz, 26 out 
of the 100 profiles). We compare the pulse profile morphologies 
with those obtained around 350 and 1400 MHz with the WSRT 
and Lovell telescopes, respectively, to study their evolution with 
respect to a magnetospheric origin and DM-induced variations. 
We do not discuss here profile evolution due to the effects of 
scattering in the ISM, which will be the target of a future work 
(Zagkouris et al. in prep.). In Sect. 2 we describe the LOFAR 
observational setup and parameters. In Sect. 3 we describe the 
analysis. In Sect. 4 we discuss the results, and in Sect. 5 we con¬ 
clude with some discussion on future extensions of this work. 


2. Observations 

The observed sample of pulsars was loosely based on a selec¬ 
tion of the brightest objects in the LOFAR-vis ible sky (declina¬ 
tion > -30°), using the ATNF Pulsar Catalog (iManchester et al.l 
l2005h for guidance. Because pulsar flux and spectral indices 
are typically measured at higher frequencies, we also based our 
selection on the previousl y published data at low frequencies 
(iMalov & Malofee"^ 1201 Oh . Since the LOFAR dipoles have a 
sensitivity that decreases as a function of the zenith angle, all 
sources were observed as close to transit as possible, and only 
the very brightest sources south of the celestial equator were ob¬ 
served. 

We observed 100 pulsars using the high-band antennas 
(HBAs) in the six central ‘Superterp’ stations (CS002-CS007) 
of the LOFAR cor^l. The observations were performed in tied- 
array mode, that is, a coherent sum of the station signals us¬ 
ing appr opriate geometrical an d instrumental phase and time de¬ 
lays (see IStanners et alJl201ll for a detailed description of L O- 
FAR’s pulsar observing modes and Ivan Haarlem et alJl201^ for 
a general description of LOFAR). The 119-167-MHz frequency 
range was observed using 240 subbands of 195 kHz each, syn¬ 
thesised at the station level, where the individual HBA tiles were 
combined to form station beams. Using the LOFAR Blue Gene/P 
correlator, each subband was further channelised into 16 chan¬ 
nels, formed into a tied-array beam. The linear polarisations 
were summed in quadrature (pseudo-Stokes I), and the signal 
intensity was written out as 245.76yus samples. The integration 
time of each observation was at least 1020 s. This was chosen to 
provide an adequate number of individual pulses, so as to aver¬ 
age out the absolute scale of the variance associated with pulse- 
phase ‘jitter’ to the cumulative profile. The jitter, also termed 
st ochastic wide-hand im pulse modulated self-noise (SWIMS, as 
in lOslowski et al.ll2()Tlll . is the variation in individual pulse in- 
tensi ty and position wit h respect to the average pulse profile (see 
also ICordesiri993l and iLiu et al.ll2012l and references therein). 
This variation does not significantly affect the measurements 
that have been carried out for the scope of this paper (i.e. pulse 
widths, peak heights), but we have checked that the resulting 
profile was stable on the considered time scales by dividing each 
observation into shorter sections and comparing the shapes of 
the resulting profiles with the overall profile. In the cases where 
stability was not achieved, we used longer integration times. Re¬ 


gardless, in almost all cases the profile evolution with observing 
frequency is a significantly stronger effect at low frequencies. 

Twenty-six of the brightest pulsars were also observed using 
the Superterp low-band antennas (LBAs) in the frequency range 
15-62 MHz. To mitigate the larger dispersive smearing of the 
profile in this band, 32 channels were synthesised for each of the 
240 subbands. The sampling time was 491.52/rs. The integration 
time of these observations was increased to at least 2220 s to 
somewhat compensate for the lower sensitivity at this frequency 
band (e.g. because of the higher sky temperature). 

For some sources, 17-minute HBA observations with the Su¬ 
perterp were insufficient to achieve acceptable signal-to-noise 
(S /N) profiles. For these, longer integration times (or more sta¬ 
tions) were needed. Hence, some of the pulsars presented here 
were later observed with 1 hr pointings as part of the LO¬ 
FAR Tied-Array All -Sky Survey for pu lsars and fast transients 
(LOTAASB see also lCoenen et al.l[2014ll . which commenced af¬ 
ter the official commissioning period, during Cycle 0 of LOFAR 
scientific operations, and is currently ongoing. LOTAAS com¬ 
bines multiple tied-array beams (219 total) per pointing to ob¬ 
serve both a survey grid as well as known pulsars within the 
primary field-of-view. 


3. Analysis 

The LOFAR HDF 50 (Hierarchical Data Format5, see e.g., 
lAlexov et alJ 1201^ data were converted to PSRFITS format 
( Hotan et al.ll2004i) for further processing. Radio frequency in¬ 
terference (RFI) was excised by removing affected time inter¬ 
vals and frequency chan nels, using the tool rfifind from the 
PRESTCO software suite (iRansotnl 200lh . The data were dedis- 
persed and folded using PRESTO and, in a first iteration, a ro¬ 
tational ephemeris from the ATNF pulsar catalogue, using the 
automatic LOFAR pulsar pipeline ‘PulP’. The number of bins 
across each profile was chosen so that each bin corresponds to 
approximately 1.5 ms. 

The profiles obtained with the HBA and, where available, 
LBA bands were compared with the profiles obtained with the 
WSRT at ~ 350MHz (from here onwards ‘P-band’) and at 
~ 1400 MHz, or with the Lovell Telescope at the Jodrell Bank 
Observatory at ~ 1500MHz (from here onwards ‘L-band’). The 
WSRT observations that we used were performe d mostly be¬ 
tween 2003 and 2 004 (see Weltev rede ^ al.ll2006l ‘WES’ from 
here onwards, and I Weltevrede et al.l2007l fOTdetailsl. The Lovell 
observations were all contemporary to LOFAR observations, 
therefore in the cases where both sets of observations were avail¬ 
able, we chose the Lovell ones because they are closer to or 
overlap the epoch of the LOFAR observations. At L-band we 
used Lovell observations for all but three pulsars: B0136H-57, 
B0450-18, and B0525H-2L In a handful of cases, where no pro¬ 
file at P-band was available from WES, we used the data from 
the European Pulsar Network (EPN) databas^B 

To attempt to align the data absolutely, we generated 
ephemerides that spanned the epochs of the observations that 
were used. This did not include those from the EPN database, 
however. Ephemerides were generated from the regular moni¬ 
toring observations made with the Lovell Telescope. The times 
of arrival were generated using data from an analogue filterbank 
(AEB) up until January 2009 and a digital filterbank (DEB) since 


* http://www.atnf.csiro.au/people/pulsar/psrcat/ 

^ The full LOFAR Core can now be used for observations and pro¬ 
vides four times the number of stations available on the Superterp (and 
a proportional increase in sensitivity). 


^ http://www.astron.nl/lotaas/ 
http://www.hdfgroup.org/HDF5/ 

^ https://github.com/scottransom/presto 
® http://www.epta.eu.org/epndb/ 
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Fig. 2. Gaussian modelling and resulting fit residuals for the pulse profile of PSR B1133+16 at four frequencies, between 0.4-0.6 in pulse phase. 
Top left; LOFAR FIBA. Bottom left: LOFAR LBA. Top right: 1.4GFIz. Bottom right: 350MFIz. The Gaussian components contributing to the fit 
are shown in colours, while in black, overlapping the profile contours, we show the best fit obtained from them. It is evident that this standard 
‘double’ profile (see Sect. ED is not well fit by only two Gaussians. Even adopting a higher number of components, the scatter o f the residuals 
is not at the same level as the off-pulse residuals (this is the criterion that was adopted for the determination of a good fit, followine lKramer et al.l 
Il994h . Nonetheless, when the residuals on-pulse were discrepant at the level of only a few percent from the ones off-pulse, we chose not to add 
more free parameters to the fit. 


then, with a typical observing cadence between 10 and 21 days. 
The observing bandwidth was 64 MHz at a central frequency of 
1402 MHz and approximately 380 MHz at a central frequency of 
1520 MHz for the AFB and DFB, respectively. The ephemerides 
were generated using a combination of PSRTIMeQ and TEMPO, 
and in the case of those pulsars demonstrating a high degree of 
timing noise, up to five spin-frequency derivatives were fit to en¬ 
sure white residuals and thus good phase alignment. 

The L-band profiles were generated from DFB observations 
by forming the sum of up to a dozen observations, aligned using 
the same ephemerides used to align the multi-frequency data. We 
re-folded both the LOFAR and the high frequency data sets using 
this ephemeris. In general, where Lovell data were available, the 
ephemeris was created using about 100 days’ worth of data. For 
the WSRT observations, an ephemeris was created spanning, in 
some cases, ten years of data and ending at the time of the LO¬ 
FAR observations. The timing procedure was the same as for 

’ http://www.jb.man.ac.uk/~pulsar/observing/progs/psrtime.html 


the shorter data spans, except that astrometric parameters were 
fit and typically more spin-frequency derivatives were required. 
The epoch of the WSRT observations is specified in Table IbTI 
in the Notes column. Given the method we used to align the pro¬ 
files, the timing solution is less accurate over these longer time 
spans than those constructed to align the Lovell data, but this too 
represents a good model, with a standard deviation (r.m.s.) of 
the timing residuals < 1 ms. We aligned the profiles in absolute 
phase by calculating the phase shift between the reference epoch 
of the observations and the reference epoch of the ephemeris and 
applying this phase shift to each data set. . 

Some of the pulsar parameters derived from these 
ephemerides are presented in Table IB. 11 The first column lists 
the observed pulsars, the second and third columns list the spin 
period and period derivative of each pulsar, the fourth column 
is the reference epoch of the rotational ephemeris that was used 
to fold the data, and Cols. 5 and 6 list the epochs of the LO¬ 
FAR HBA and LBA observations. In Cols. 7 and 8 two mea- 
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Table 1. Pulsars for which the absolute alignment was not achieved 
with the refolding using the same ephemeris (see text for details). The 
extra DM shift (in cm^^pc) and corresponding phase shift (A0) needed 
to align the profiles are indicated, or possibly other reasons for the ob¬ 
served shift, e.g., S for scattering, which notably alters the shape of the 
profile, g in case a glitch occurred during the range of the ephemeris, 
number of spin-frequency derivatives fitted to obtain a good ephemeris 
(i.e. F#), or final r.m.s. (in ms) of the best timing solution. 


PSR Name 

extra DM shift/causes for misalignment 

B0114-H58 

S 

B0525-H21 

F4, S, g 

B1633-H24 

F4, rms=1.3, DM=-0.11, A0=O.O43 

B1818-04 

S 

B1839-H09 

rms=1.4, DM=-0.13, A0 = 0.051 

B1848-H13 

DM=-0.04, A(P = 0.017 

B1907-H10 

F3, S 

B1915-H13 

S 

B2148-H63 

S 


surements for the DM are given: the first as originally used to 
dedisperse the observations at higher frequencies, and the sec¬ 
ond as the best DM obtained from the fit of the HBA LO FAR 
observations using PRESTO’S prepfold (iRanso ml l200l . The 
next three columns provide the pulsar’s spin-down age, mag¬ 
netic field strength, and spin-down luminosity as derived from 
the rotat ional parameters accordin g to standard approximations 
(see e.g. lLorimer & Krameiil200^ : 

t[s] = 0.5P/P, (2) 

B[G] = 3.219-10'“’V^, (3) 

£[erg/s] = • 10^^PlP\ (4) 

where P is measured in s and P in ss '. The resulting aligned 
profiles for the 100 pulsars are shown in Fig. IBT] in the ap¬ 
pendix. The star next to the name of the band (P-band in most 
cases, with the exception of B0136-I-57, where we used the P- 
band for absolute alignment) indicates that the corresponding 
band was aligned by eye based on the absolute alignment be¬ 
tween the LOFAR data and the other high-frequency band. The 
alignment was made based only on the choice of a specific point 
along the rotational phase of the pulsar, at the reference epoch of 
the ephemeris, but unmodelled DM variations can be responsible 
for extra, albeit small, phase shifts (up to a few percent, see Ta¬ 
ble [T]i. Indeed, observations performed at different times, quite 
far apart, and at different frequencies, can possess quite differ¬ 
ent apparent DMs (up to some tenth of a percent, see Table [U. 
DMs that are due to the ISM have, as expec ted, a time depen¬ 
dence (see I You et alll2007i iKeith et ^l2()13h . and these differ¬ 
ences can become quite relevant especially at the lowest LOFAR 
frequencies. We chose to re-dedisperse all the profiles, LOFAR 
and high-frequency ones, using the DM obtained as the best DM 
with prepfold for the HBA LOFAR observations, prepfold 
determines an optimum DM by sliding frequency subbands with 
respect to each other to maximise the S /N of the cumulative pro¬ 
file. The intra-channel smearing caused by DM over the band¬ 
width at the centre frequency is indicated by the filled rectangle 
next to each profile in Fig. IBT] 


Table 2. Pulsars for which one or more glitches occurred during the 
range of the ephemerides used in this paper (above the horizontal line) 
or near the range of validity of our ephemerides (below the horizontal 
line). The uncertainty of the glitch epoch is quoted in parentheses and it 
is expressed in days, except for the case of PSR B1907-t-10, where it is 
quoted in seconds. Sources: Jodrell Bank and ATNF glitch archives. 

References. 1 lljEs pinoza et alj ll201 Ih : 12l lJanssen & Stannersl l l2006h : 
131lYuanetal] 112010). _ 


PSR Name 

Start 

[MJD] 

End 

[MJD] 

Glitch Epoch 
[MJD] 

B0355-)-54 

51364.6 

56262.2 

51679(15)['l 

51969(1)[" 

52943(3)['l 

53209(2)['l 

B0525H-21 

52274.9 

56641.1 

52296(1)^-3' 

53379PI 

53980(12)P1 

B0919H-06 

55555.2 

56557.5 

55152(6)'*' 

B1530H-27 

51607.0 

56535.6 

49732(3)"' 

B1822-09 

54876.5 

56571.8 

54114.96(3)"' 

B1907H-00 

54984.1 

56556.9 

53546(2)"' 

B1907H-10 

54924.5 

56535.1 

54050(350 [s])[3] 

B2224H-65 

55359.7 

56570.1 

54266(14)[2-3] 

B2334H-61 

54635.4 

56507.3 

53642(13)"'3' 


Only in a few cases, documented in Table [T] was the DM 
obtained from the LOFAR observations not used for the align¬ 
ment. Those are the cases for which, also evident from Fig. IBTI 
the intra-channel smearing caused by DM is a substantial frac¬ 
tion of the profile width (similar to or higher than the on-pulse 
region), and therefore the quality of the measurement is lower 
than that obtained at a higher observing frequency. On the other 
hand, in some other cases (although rare, see Table IB. Il l, even 
the LBA measurement was good enough to provide a DM mea¬ 
surement, and in these cases we were able to use that for the 
alignment. In this way we obtain the best alignments, in general, 
even though some residual offsets could still be observed in a 
handful of cases. 

For those cases (listed in Table [T]i where the remaining off¬ 
set was noticeable by eye, we investigated the possible causes 
after refolding and applying the new DM. We checked whether 
the pulsars in our sample had undergone any glitch activity dur¬ 
ing the time spanned by our ephemerides. Sixteen out of our 
100 pulsars have shown glitch activity at some time. Seven of 
them have experienced glitches relatively close to the epoch of 
our observations, but only two of them during the time spanned 
by our ephemerides. The relevant glitch epochs of these pulsars 
are presented i n Table El and were taken fr om the Jodrell Bank 
glitch archival (lEsnino^ et al.ll2()Tll l2012l) . integrated with the 
ATNF pulsar archiv^B- We note that while the glitch activity 
could have had an influence on the shift of PSR B0525-I-21, the 
recurrent activity of PSR B0355-1-54 did not cause as notable an 
impact on the alignment. In some cases the profile is scattered in 
the LOFAR HBA band and rapidly becomes more scattered to¬ 
wards lower frequencies. This will affect the accuracy of the DM 
measurement, potentially causing an extra profile shift between 

® http://www.jb.man.ac.uk/pulsar/glitches/gTable.html 
^ http://www.atnf.csiro.au/people/pulsar/psrcat/glitchTbl.html 
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bands. In some cases the ephemeris had a large r.m.s. timing 
residual, or it included higher order spin-frequency derivatives 
beyond the second, which is typical of unstable, ‘noisy’ pulsars. 
All these cases are indicated in Table [T] In the other cases, we 
calculated that a DM difference < 0.05 cm^^pc, compatible with 
the measurement uncertainties, would compensate for that shift. 
For this reason, we applied an extra shift to align these pulsars’ 
profiles by eye in Fig. IB. H and referenced this shift in Table[T] 
Only a small number of pulsars in our sample show in¬ 
terpulses: B0823-H26, B0950-H08, B1822-09, B1929 h- 10, and 
B2022-I-50. For these pulsars the profile longitude is shown in 
its entire phase range instead of in the phase interval 0.25 - 0.75, 
as was chosen for the other profiles. The phase-aligned pro¬ 
files, which normally have their reference point at 0.5, have been 
shifted by h- 0.25 in these cases, to help the visualisation of the in¬ 
terpulse. A zoom-in of the interpulse region is shown in Fig. IRTI 
(LBA data were excluded because none of the interpulses were 
detected in that band). Our sample also contains a few moding 
pulsars (notably B0823-I-26 and B0943-I-10). For these pulsars 
we caution that the profile reflects only the mode observed in 
the particular observations presented here. A more detailed treat¬ 
ment of the low -frequency pr ofiles of B0823H-26 and B0943H-10 
can be found in ISobev et al.l (1201 5l) an d iBilous et al.l (l2014l) . re¬ 
spectively. In yet other cases, for instance B1857-26, no moding 
behaviour has previously been documented, but notable changes 
in the profile are seen from high to low frequencies. These might 
reflect different modes of emission when the various observa¬ 
tions were taken. A more detailed discussion of some cases of 
peculiar profile evolution can be found in Sect. 15.31 

We fit the multi-band profiles of each pulsar using Gaussian 
components (see an example in Fig.|2]i, whic h are a good repre- 
sentation for the profiles of slow pulsars ('see lKramer et al.ll 19941 
and references therein). We used the program pygaussfit. py, 
of the PRESTO suite, which has the advantage of providing an 
interactive basis for the input parameters to the Gaussian fits. 
This program can be used to apply the same method as in 
iKramer et al.l (1 19941) (see e.g. Fig. 3 of their paper), as it provides 
post-fit residuals (for the discrepancy between the model and the 
data, see bottom half-plots of Fig. |2]i that we required to have 
approximately the same distribution in the on-pulse as in the off- 
pulse region. The full rotational period, and not only the on-pulse 
region, was taken into account by the fit, also allowing for distin¬ 
guishing interpulses or small peaks at different phases from the 
noise. The Gaussian components derived using this method were 
chosen to satisfy the condition of best fit with minimal redun- 
dancy, and no physi cal significance should be attributed to them. 
iHassall et al.l (l2012h have shown that it is possible to accurately 
model the evolution of the profile with frequency using Gaus¬ 
sian components, but a specific model has to be applied individ¬ 
ually to each pulsar, requiring much careful consideration. Such 
an analysis is beyond the scope of this paper. In the absence of a 
more comprehensive treatment of scattering for LOFAR profiles, 
which is envisaged for subsequent papers, the most evident scat¬ 
tering tails of the low-frequency profiles have not been modelled 
and the corresponding profiles components were not included in 
the table and are not used in any further analysis. 

We used the mathematical description of the profiles in terms 
of the Gaussian components to calculate the widths and ampli¬ 
tudes of the observed peaks. For each profile we obtained the 
full width at half maximum (wso) and the full width at 10% of 
the maximum, wio- To calculate the errors on the widths, we 
simulated 1000 realisations of each profile, using the noiseless 
Gaussian-based template and adding noise with a standard de¬ 
viation equal to that measured in the off-pulse region of the ob¬ 


served profile. By measuring the widths in these realisations, we 
obtained a distribution of allowed widths from which we deter¬ 
mined their standard deviation. We note that these errors are sta¬ 
tistical only and do not take into account the validity of our as¬ 
sumptions, that is, the reliability of the template used. To cross¬ 
check the width of the full profile, we tried different methods. 
An example of how these widths were calculated is described in 
Sect.l^of the Appendix and is shown in Fig. IA.il In Table IB. 2l 
Cols. 2 and 3 show wso and wio in degrees at all frequencies, 
calculated using the Gaussian profiles and cross-checked using 
the on-pulse visual inspection. The last two columns represent 
the calculated spectral index b of the evolution of these widths 
(wso and wio) with frequency as ic oc v^. The table shows that this 
fit can be highly uncertain. We also note that for more data, the 
linear fit is not always the best representation of the real trend of 
the profile evolution (see Fig. |3] and a more detailed discussion 
in Sect. l5.1.TT i. Column 4 of Table lB.2l lists the duty cycle of each 
pulsar, calculated as wwlP, where P is the pulsar period. 


4. Results 

Here we present the results of LOFAR observations of 100 pul¬ 
sars, considering their profile evolution with frequency and com¬ 
paring them with observations at higher frequencies. In particu¬ 
lar, we study how the number of profile peaks, their widths, and 
the relative pulse phases vary with frequency. 


4.1. Pulse widths 

We calculated the evolution of the profile width across the fre¬ 
quency range covered by our observations. We chose icio for our 
calculations and cross-checked using Wop, as it is better suited 
for multi-peaked profiles than and less affected by low S fN 
than Wpow (see Sect. 0 Fig. IA.il and Table IB. 211 . When the mea¬ 
surements disagreed, the Gaussian fit was refined after visual in¬ 
spection of the obtained width. 

We calculated the dependence of the width of the profiles on 
the pulse period, considering the different frequencies separately. 
We note that the pulse width is not a direct reflection of the beam 
size or diameter (i.e. 2p, where p is beam radi us). For a visual 
repres enta tion of the geometry see for instance iMaciesiak et al.l 
(l20TTli and iBilous et al.l(l2014h . In fact, only if the observer’s line 
of sight cuts the emission centrally for magnetic inclination an¬ 
gles, O', that are not too small (i.e. a >~ 60°), w » 2p. In such a 
case, when the emission beam is confined by dipolar open held 
lines, we would expect a dependence, which has indeed 
been observed when correcting for geometrical effects by trans- 
forming the pulse width into a beam radius measurement (see 

lRankirill993l:lGil & Krawczvkll996l:lMaciesiak et al.l2012h . For 

circular beams, profile width and beam radius are related by the 
relation first derived bv iGil et^(ll984l) : 


pio = 2 sin 


sin a sin (a + ff) sin^ 



(5) 


The angle p is the impact angle, measured at the fiducial phase, 
(p, which describes the closest approach of our line of sight to 
the magnetic axis. This equation is derived under the assumption 
that the beam is symmetric relative to the fiducial phase. Typi¬ 
cally, widths are measured at a certain intensity level (e.g. 50% 
or 10%, as here), and p values are derived accordingly. In many 
cases, profiles are indeed often asymmetric relative to the chosen 
midpoint, or become so as they evolve with frequency. We note 
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Fig. 3. Width at 10% of the maximum (loio) encompassing the outer components of the profile, where present. The plot shows the evolution of 
wiQ as a function of observing frequency for each pulsar. 


Article number, page 7 ofl36l 


































































































































































































































A&A proofs: manuscript no. 25196_ap 


150 MHz 



Period (s) 


150 MHz 



Fig. 4. HBA profile widths (top) and icio (bottom) as a function of spin period. Left side: The blue and red dots represent the data (red: 
interpulse pulsars, blue: other pulsars), while the red solid line represents the best fit (non-weighted ) an d the red dashed lines represent its Icr 
dispersion for lu = A • , The fit is calculated using our interpulse pulsars, following iRankinI il990h and iMaciesiak et ahl ll201 ih (see discussion 

in the text). Right side: The blue dots represent the data, while the green solid line represents the best (non-weighted) fit for all the pulsars: 
w = A ■ P^" , and the green dashed lines represent its Icr dispersion (see discussion in the text). 


that for a central cut of the beam (ff - 0) and for an orthogonal 
rotator (a ~ 90°) the equation reduces to p = uj/2 as expected, 
while in a more general case, where ff — Q and a :s> p the re¬ 
lation reduces to p = (w/2) sin a. In principle, it is possible to 
determine a and (3 with polarisation measurements. However, in 
reality the duty cycl e of the pulse is often too s mall to obtain re¬ 
liable estimates (see iLorimer & Krameill2004l). A lternatively, at 
least for a, the relation reported bv iRankinI (Il993h can be used: 

u^50,core(lGHz) = 2.45° ■ sin(a), (6) 


calculated from the observed width dependency on period for 
the core components of pulsars (see Sect. EB, which is intrinsi¬ 
cally related to the polar cap geometry. Equation |6] is valid at 
1 GHz, but can be applied at LOFAR frequencies, maintain¬ 
ing the same dependence, if the impact angle (i <sc Pm rp.; sincr 
should be ignored for or th ogonal rotator s . Add it ionally, Ran^ 
1993h: lr.il et al.l (Il993h : iKramer et al.l (11994 : iGoiild & Evne 


19981) suggested that ‘parallel’ p — P relations are found if the 


radio emission of the pulsar can be classified and separated into 
emission from ‘inner’ and ‘outer’ cones, which seem to show 
different spectral properties (see Sect. l5.1l for details). 

Figure |4] represents the 50% and 10% widths of the pro¬ 
files as a function of the pulsar period. We show the re¬ 
sults for FOFAR data, using only the HBA data, for which 
we have the largest sample. In the l eft pa nel of Fig. 01 we 
adopted the assumption s from iRankinI (Il990h . later followed by 
iMaciesiak et alJ (1201 Ih . and used our interpulse pulsars (over¬ 
lapping their samples of ‘core-single’ pulsars that show inter¬ 


pulses) as our orthogonal rotators to calculate a minimum esti¬ 
mated width using a fixed dependence on the period as P 
The interpulse pulsars in our sample are plotted in red in Fig.0] 
and are shown in Fig. lB.TI and labelled IP in Table IB. 2l The red 
solid and dashed lines represent the best fit of the dependency of 
(Uso and wio on which should constitute a lower limit to 

the distribution of pulse widths. We obtain 

W5o(150MHz) = (3.5 + 0.6)° • p-0-^o±o.o2 
wio(150MHz) = (10 + 4)° ■ p-o so±o.o2^ 

where P is in seconds, the error is quoted at Icr for the amplitude, 
and the er ror on the power-law in dex -0.50 was taken to be 0.02, 
following lMaciesiak et al.l (1201 ih . 

In the right panel of Fig. 0] we present a fit to the widths of our 
FOFAR sample as a function of pulse period. Because the scatter 
is much larger than the individual error bars, we performed a 
non-weighted fit. The lines represent the best fit to the data (solid 
line) and its icr dispersion (dashed lines). Here we calculated 

u;5o(150MHz) = (6.2 + 0.8)° • p-o i±0-4 (9) 

i(;io(150MHz) = (16 + 2)° • T--o.3±o.4^ 
where the errors are 1 cr. 

The widths follow an inverse dependency with the pul¬ 
sar period, co nsistent with p revious analyse s at h igher fre¬ 
quencies (e.g. iRankinI 1 19901: IMaciesiak et al.l 1201 ll but also 
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150 MHz 



Period (s) 


1500 MHz 



Period (s) 

Fig. 5. Duty cycle (w[q/P) vs period of the pulse longitude at the two 
frequencies for which we have data for all the pulsars in our sample: 
HBA band from LOFAR and L-band , for comparison. The error bars 
are omitted, as in lLorimer et al.l (Il99.‘5h . to more easily guide the eye on 
the trend. 


iGil et al. 1^93; Arzoumanian et al1l2002l) and at these frequen¬ 
cies ( Kuzmin & Losovskvlll999(r in general, broadening by ex¬ 
ternal effects may also be expected, even though not dominant: 
while we were careful to avoid evidently scattered profiles in our 
sample, DM smearing can also contaminate it. Finally, because 
our data are chosen according to detectability of the pulsars at 
low frequencies, a different bias in the observed sample com¬ 
pared with high frequencies cannot be excluded. In conclusion, 
our determination of the relationship between W 50 or wiq and 
P is only a first step to determining the relation for the model- 
independent beam shape, which is to be determined when more 
polarisation measurements are available. 


Figure |5] presents the duty cycle (wio/F) of the pulsars in 
two bands: LOFAR HBA and L-band, for comparison, plotted 
against the period. The values of the duty cycle are reported in 


Col. 5 of Table IB. 2I The inverse correlation that is observed, im¬ 
plying that shorter period pulsars have larger beams, is also evi¬ 
dent at LOFAR frequencies. It can be used to characterise pulsar 
beams and help create accurate beaming models for pulsars in 
the Milky Way, which wo uld in turn constrain the Galactic pop¬ 
ulation and its birth rates (iLorimer et al.lll99'5h . 

4.2. Spectral evolution of individual components 

No absolute flux calibration of beam-formed LOFAR data was 
possible with the observing setup used for these observations. 
Therefore no spectral characterisation could be attempted yet. 
Nonetheless, we attempted a characterisation of the relative am¬ 
plitudes of pulse profile components for pulsars with multiple 
peaks. In Table IB.3l we list the pulsars for which double or mul¬ 
tiple components can be observed and separated in at least two 
frequency bands. For these pulsars we selected the two most 
prominent peaks and calculated the evolution of the relative 
heights with observing frequency. We chose to select the peaks 
as the two most prominent maxima of the smoothed Gaussian 
profiles and verified by eye that we were consistently follow¬ 
ing the same peak at all frequencie s. We note that as a result 
of subtle profile evolution (see e.g. iHassall et alJl2012h . shifts 
of the peaks in profile longitude cannot be excluded, which we 
did not track here. The profile evolution is in some cases quite 
complex and the profiles are sometimes noisy, therefore the pro¬ 
files as presented in Fig. IBT] need to be reviewed before draw¬ 
ing any strong conclusions based on Fig. |6] We first attempted 
to apply a power-law fit to calculate how the ratio between the 
two peak amplitudes changes with observing frequency for all 
the pulsars: [P2/Pl]{v) oc v^, where PI is the peak at the ear¬ 
lier phase and P2 is the peak at the later phase. This fit pre¬ 
sented large errors, and the distribution of the spectral indices 
was peaked close to 0, with an average of 0 + 2, indicating no 
systematic evolut ion despite the larg e scatter. A similar finding 
was obtained by IWang et aH (l2001l) . who only selected conal 
double pulsars. They concluded that a steeper spectral index for 
the leading or trailing component are equally likely, arguing in 
favour of a same origin of the peaks in the magnetosphere, as 
expected if both components correspond to two sides of a conal 
beam. They also found a dominance of small spectral indices, 
with a quasi-Gaussian distribution, indicating no systematic evo¬ 
lutive trends. The different evolution of the peaks would then be 
due to geometric beaming effects. While in their case the pulsars 
were carefully selected so as to include only the conal doubles, 
in our case no such distinction was followed, so that the different 
relative spectral indices could also depend on a different origin 
of the emission regions (see also Sect. I53l l. 

Given the large scatter of this result, and because our sample 
is heterogeneous (with double and multi-peaked pulsars), we in¬ 
vestigated the ratios more closely. Figure |6] shows the evolution 
of the ratios with observing frequency for each of the pulsars 
used for this calculation. We ordered the pulsars into two groups, 
taking first the pulsars for which the profiles were already stud¬ 
ied in previous works, and sorted by right ascension within the 
groups. In most cases it is apparent that the simple power law is 
not a good fit to the data and can be misleading if measurements 
are only possible at two frequencies. 

Table IB.3I provides a detailed summary of these measure¬ 
ments. We observe that the prominence of the peaks seems to 
shift from low to high frequencies with, in most cases, a net in¬ 
version in the dominance of FI from LOFAR HBAs to L-band. 
The inversion point is also indicated in Fig.^by the blue hori¬ 
zontal line. In general, starting from LOFAR frequencies, there 
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Fig. 6. Ratio of peak amplitudes for pulsars with multiple peaks as a function of observing frequency. The red line connects the dots (it is not 
representative of the power law used to fit an exponential dependence: P2/Pl(v) oc It is evident that in most cases a power law does not represent 
the best fit for the frequency evolution of the peaks’ ratio, although care should be taken as it is hard to reliably track the peak amplitudes in some 
cases and to follow the same P2 and PI (see Sect. lS.l.TI for details). The blue line corresponds to the inversion point where the peaks have equal 
amplitude. It is evident that in a number of cases the relative amplitudes of the peaks invert as a function of frequency (see Sect. lS.l.TI for details). 
The black vertical lines in between plots 11 and 12 mark the change between previously studied cases and new (see text for details). 


seems to be a trend that the peaks change from being more dis¬ 
similar amplitudes at low frequencies to bec oming more s imilar 
at high frequencies. We note that Fig. 1 of IWang et al.l (l200lh 
shows that a linear trend of the peaks’ ratio with frequency fits 
the data well in most cases, meaning that at higher frequen¬ 
cies the relative amplitude of the peaks will again depart from 
equality. Notable changes in the observed pulse profile prop¬ 
erties at low frequencies with respect to high frequencies hav e 
previously been observed for i nstance by iHa ^saUet^ (l2012h . 
iHankins & Rankinl(l2010l) . and lizyekoya et alJ (Il993h . 

An obser yed feature that can contribute to this behayiour was 
discussed by lHassall et al.l (l2012h : they modelled the profile eyo- 
lution with Gaussian components that were free to eyolye longi¬ 
tudinally in a dynamic template. The examples presented there 
(two of which are also in our sample: B0329+54 and B1133+16) 
showed that the components change amplitudes and moye rela- 
tiye to one another. 

Hardening of the spectrum of the second peak is obseryed in 
gamma-rays in th e typical case of two prominent caustic peaks 
(I Abdo et ^l20 1 3h and is explained with the different paths that 


curyature-emitting ph otons follow in t he leading and trailing side 
of the profile (see e.g. lHirotanil|201 Ih . but it is not obyious that 
this should also follow for the radio emission. 

5. Discussion 

5.1. Phenomenological models for radio emission 

Based on the findings discussed aboye, we drew some conclu¬ 
sions on the models that haye been proposed to explain the ob¬ 
seryed properties of pulsar profiles and on some predicted ef¬ 
fects such as radius-to-frequency mapping (RFM). These models 
haye largely been developed based on observations performed at 


> 200 MHz. 


Rankin’s model 






Radhakrishnan & Rankinll 19901 

RankinI 19931 iMitra & RankinI 


2002 h proposed that the emission comes from the field lines 


originating at the polar caps of the pulsar, forming two con¬ 
centric hollow emission cones and a central, filled, core. There 
is a one-to-one relation between the emission height and the 
observing frequency, so that at different frequencies the profile 
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evolves, as more components come into view or disappear. 
Rankin based her classification on the number of peaks and 
polarisation of the pulse prohles. P rohles with u p to hve 
components are observed (although see iKrami^l 1 9941 for a more 
detailed classihcation). The prohles can be single (5, or Sd 
based on whether the prohle will become triple or double at 
higher or lower frequencies, respectively), double (D), triple (T 
or cT), tentative quadruple (cQ), and quintuple (indicated as 
multiple M), where c represents the presumed core origin. 

iLvne & Manchester (Il988h found that their data agreed with 
the hollow cone model, and they also observed a distinction in 
spectral properties between core and cone emission, or at least 
inner and outer emission. However, based on asymmetries of 
the components relative to the midpoint of the prohle and the 
presence of so-called partial-cone prohles, they proposed that a 
window function dehnes the prohle shape, and within this, the 
locations of emission components can be randomly distributed. 

A further step i n this direction was made by 
iKarastergiou & .lohnstonl (l2007h . who assumed a single 
hollow cone structure without core emission but instead with 
patches of emission from the cone rims. Emission could come 
from different heights at the same observing frequency, but still 
following REM; the number of patches changes with the age of 
the pulsar: up to ~ 10 patches, but only at one (large) height for 
the young pulsars, and up to ~ 4 - 5 at ~ 4 different (lower) 
heights for the older pulsars. This also explains the narrowing 
of the prohles as the period increases (see Sect. 14.11) . and their 
simulations successfully reproduce the observed number of 
prohle components (i.e. typically Ncomp ^ 5). The central 
components are then simply more internal and surrounded by 
the external ones coming from higher up in the magnetosphere: 
this is why they show single peak prohles and steeper spectra. 
Younger pulsars have been observed to have simpler prohles, 
but typically with longer duty c ycles than those of older ones. 
IKarastergiou & JohnstonI (l2007h predicted that there should be 
a maximum height of emission of ~ 1000 km. The minimum 
height, on the other hand, is quite varied but is close to the 
maximum allowed for young pulsars, which then emit only from 
one or two patches. Because of the width depends on the period, 
the opening angle of the cone would then be comparatively 
larger at the same height for younger than for older pulsars. 

5.1.1. Radius-to-frequency mapping 

In the framework of the standard models for pulsar emission, 
where the radio emission is predicted to come fro m the polar 
caps of t he pulsar, it has been hypothe sised (e.g. iKomesar^ 
Il970l a nd[Ruderman & Sutherlandlll975h and in some cases ob¬ 
served (ICordeslll978l) that the emission cone widens as we ob¬ 
serve it at lower frequencies because we are probing regions fur¬ 
ther away from the stellar surface where the opening angle of the 
closed magnetic held lines is broader. The phenomenon is more 
apparent at low frequencies (< 200MHz) and is therefore ideal 
to study using LOEAR. 

A limited observational sample h as always biased t he con¬ 
clusions about REM. Originally (e.g. lKomesarofflll970t) it was 
thought that the REM behaviour could be observed as a power- 
law dependence of the increase in peak separation with decreas¬ 
ing observing frequency, and asymptotically approaching a con¬ 
stant separation at high frequencies (> 1 GHz). It was therefore 
proposed that two power laws (i.e. two different mechanisms) 
regulated the evolution of the pulse prohle, with a break fre¬ 
quency at approximately 1 GHz. 


la orsetd (Il99ll) analysed this dependency using a sample of 
pulsars observed at various frequencies. He concluded that no 
break frequency seemed to be necessary to model the evolution 
of the pulse prohle. On the other hand, a simple power law (or 
a quadratic dependence, indiscernible with his data) and the ad¬ 
ditional constraint of a minimal emission width (or peak sepa¬ 
ration) could ht the data at all frequencies. He obtained the fol¬ 
lowing functional dependence from a phenomenological model 
of component separation: 

Ae ^ A ■ ( 11 ) 

where A6 is the component separation, 6 is the separation power- 
law index of the components, and the constant value at 

high frequency that the pulse separation tends to. The predictions 
for 6 are quite varied depend ing on the theoretical model (see 
Ta ble 1 inIXilouri s et al.lli 996l) . 

iMitra & RankinI (l2002ll also did extensive work on REM. 
They assumed double prohles to derive from conal emission 
and therefore analysed a sample of ten bright pulsars showing 
prominent cone components. They found that inner cones are 
not affected by REM and that their component separation does 
not vary with observing frequency, while outer cone components 
show REM and increase their separation with decreasing observ¬ 
ing frequency;_^ 

iHassall et alJ (12012h have discussed that REM does not seem 
to be at play for some pulsars observed with LOEAR. Here, 
with a more statistically signihc ant samp l e, we can discuss the 
m atter in more deta i l. Eoll owing lThorsetd (Il99lh and in particu- 
lar lMitra & RankinI (l2002h . we investigated double-peaked pulse 
prohles and their component separation (see peak phases in Cols. 
5 a ndbofTablelR^. 

iMitra & RankinI (l2002h divided a group of double prohle pul¬ 
sars into three groups The pulsars from groups A and B are asso¬ 
ciated with outer cone emission (the difference between the two 
being a ht with or without the constraint po = Ppc, where po is the 
constant equivalent to A0„,/„ from Eq.[TT] relative to the beam ra¬ 
dius, and Ppc is the beam radius at the polar cap edge) while the 
pulsars from group C are associated with inn er cone emission. 
Of the ten pulsars of IMitra & RankinI (l2002h . the pulsars from 
our sample that fall inside each group are 
Group A: B0301H-19, B0525-H21, and B1237-H25 
Group B: B0329+54 and B1133-1-16 
Group C: B0834+06, B 1604-00, and B1919H-21. 

IMitra & RankinI (l2002l) reported that pulsars from groups A and 
B show REM, while the pulsars from group C show almost 
no evolution at all. Although the prohle in our sample evolves 
rapidly at low frequency, it seems that a similar behaviour can 
be observed (see single cases in Eig. E). 

Eigures |3] and |7] show a similar calculation using tuio. We 
plotted the evolution of wio as a function of observing frequency 
for the single-peaked pulsars and the histogram of the spectral 
indices of this evolution. We excluded the LBA and HBA pro¬ 
hles that showed signihcant scattering. The errors were calcu¬ 
lated from the Gaussian ht, taking into account both the noise 
contribution and any unaccounted scattering of the prohle. Al¬ 
though the values in Eig. |3 seem to follow the power-law, the 
prohle width in some cases effectively behaves in a non-linear 
way, as can be cross checked in Eig. lB.ll for the single cases. 

The weighted mean spectral index from Eig. |7]is d = -0.1(2). 
Our results are compatib le at Icr with the predictions made by 
iBarnard & AronsI (Il986l) : the component separation does not 
vary (no REM, 6 ~ 0.0), but the distribution in Eig. |7] peaks at 
negative spectral indices, which is evidence for a weak widening 
of the prohle at low frequencies. Eollowing the predictions of 
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Fig. 7. Distribution of the indices 6 of the power law used to compute 
the evolution of Ww across our frequency range as u;io(v) oc y*. 

iGil & KrawczvkI (Il996h . the dependence of ww with observing 
frequency based on their calculations should be (5 = -0.21 for 
RFM and conal beams. This is in support of their model, while 
the fact that we see a broad distribution and a flatter median in¬ 
dex might be explained by a subdivision of pulsar behaviours 
according to the Rankin groups. For a future complete analysis, 
geometry should be taken into account to perform a study on the 
beam radii (p) rather than the pulse widths. 

An alternative but complementary explanation to the ob¬ 
served widening of pulsar profiles with decreasing observing fre¬ 
quency can be found in the theory of birefringence o f two differ¬ 
ent pr opagation modes of a magneto-active plasma (iMcKinnonI 
Il997h . These two modes of propagation follow a different path 
along the open held lines. The nature of birefringence is such 
that the two polarisations are spatially closer together at higher 
frequencies, and de polarisation will occur where they overlap. 
iBeskin et al.l (Il988l) predicted that the two modes of propagation 
would result in two different indices: 5 = -0.14 and 6 = -0.29 
for the ordinary mode and d = -0.5 for the extraordinary 
mode. While our observations are compatible with either sce¬ 
nario or a combination of them, polarisation studies will help 
discern between the two interpretations of this phenomenon (see 
iNoutsos et al.ll2015l) . 

5.1.2. Profile complexity 

iKarastergiou & JohnstonI (l2007h searched for a relation between 
the number of peaks in the profile of a pulsar and some observed 
or derived parameters, such as its period, age, and rotational en¬ 
ergy loss. As a general trend, they found that faster, younger, 
more energetic pulsars would typically show less complex pro¬ 
files, which prompted them to assume that the regions of emis¬ 
sion for these pulsars arise at higher altitudes in the magneto¬ 
sphere and are, therefore, less numerous (see Sect. |5TT) . Notably, 
an abrupt change in this respect can be observed at P < 150 ms, 
T < 10^ yr and E > 10^^ ergs '. 

We calculated the same relations using our sample of 100 
pulsars and made the comparisons using LOFAR HBA band and 
the L-band data. Figure [8] shows the relation between the pe¬ 
riod of the pulsar and its spin-down age r for the two frequency 
bands. Each circle represents a pulsar, and its colour and diam¬ 
eter represent a different number of peaks of its profile. The cir¬ 


cles are larger for growing number of peaks with frequency, and 
the number of peaks is ordered by colour, in the order white, 
red, green, blue, and purple. In the histograms we summed the 
pulsars separated by number of peaks in their profiles to search 
for trends as a function of either period or t. No trends are ev¬ 
ident in any of the histograms. Our sample does not cover the 
region of young energetic pulsars in a statistically significant 
way so that while the few cases might confirm the predictions of 
IKarastergiou & JohnstonI (l2007l) . nothing in favour or disfavour 
can be stated in this respect. 

Comparing the histograms from the two plots it appears, in 
general, that the number of peaks changes from one band to the 
other. For more clarity, we plotted the trend of the number of 
peaks from HBAs to L-band in Fig.|9] Here the triangles pointing 
down (red) represent the pulsars that show fewer peaks at HBA 
frequencies than at L-band: a total of 30 pulsars; and the ones 
pointing up (blue) are those for which the number of peaks is 
higher at HBA frequencies than at L-band: 16 pulsars. The green 
triangles pointing right represent the pulsars whose number of 
peaks does not change between the two frequencies: 30 pulsars. 
There is a slight indication that faster-spinning pulsars have a 
higher number of peaks at L-band than at HBA frequencies while 
the slower spinners have more peaks at HBA frequencies than at 
L-band. The younger pulsars in general have more peaks at L- 
band than HBAs; the older have fewer peaks at L-band than at 
HBAs. In both cases the null hypothesis probability that the two 
distributions are the same based on the Kolmogorov-Smirnoff 
statistics cann ot be rejected at the 20% con fidence level. This, as 
pointed out bv iLvne & Manchesteil(ll988h . might be a selection 
effect or a resolution effect (as fast pulsars tend to have wider 
pulses that can be resolved more easily), but given our statistics, 
we can also tentatively assume that our sample includes pulsars 
with different behaviours. 

Additionally, a number of pulsars show an increase in the 
number of components at HBA compared to L-band, which may 
also be explained at HBA frequencies by the general expectation 
that a higher portion of the beam can come into view at lower 
frequencies, according to RFM. On the other hand, the fastest 
and youngest pulsars show an opposite trend , similar to what 
is als o observed in millisecond pulsars (see iKondratiev et al.l 
12015h . which are even faster and therefore have a wider duty- 
cycle (see Fig. |5] l. A si milar finding was also reported by 
iHankins & Ricketfl (Il986l) . who analysed the frequency depen¬ 
dence of pulsar profiles for 12 pulsars in the 135-2380 MHz fre¬ 
quency range. They argued that the occurrence of single profiles 
at low frequencies that become multiple at high frequencies can 
be explained in the framework of the ‘core and cones’ models (in 
the formulation of lRankinlll993ll . In particular, this is expected 
to occur if we can only observe the core emission at low fre¬ 
quencies, which is observed to typically have a steeper spectrum, 
while the outriding con al components only e merge at higher ob¬ 
serving frequency (see iKramer et al.l 1 19941 for arguments why 
this is caused by geometrical reasons and applies to inner and 
outer components regardless of their ‘nature’ as core or cones). 
A word of caution is needed here, related to how the number of 
peaks was determined: it is possible that we achieved a good fit 
to the profile using a smaller number of Gaussians in HBAs rel¬ 
ative to L-band because the quality of the profile is lower and 
so fewer components need be fit (for details on the method see 
Sect. [2 and the single cases can be studied by comparing Lig. 
iBTland TableESi. 
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Fig. 8. Spin-down age vs period in FlBAs and at L-band for the pulsars 
whose profiles are not affected by scattering in HBAs. The colour (and 
symbol size) code for the number of peaks is white=l, red=2, green=3, 
blue=4, and violet > 4. The histograms at the sides represent the distri¬ 
bution of these values as a function of spin-down age (right-most his¬ 
togram) and of period (top-most histogram). From the plots no trend is 
visible in the number of peaks. 


5.2. External effects on profile evolution 

The interstellar medium affects the pulse signal while it trav¬ 
els towards the observer, and it strongly depends on observ¬ 
ing frequency, with observations at low frequencies being more 
strongly affected by sc attering and dispersion delay (see also 
IZakharenko et al.l2013h . We here did not correct the profiles pre¬ 
sented for scattering effects that can smear out the signal espe¬ 
cially at LBA frequencies, but we considered the effect of intra¬ 
channel dispersive smearing. 



Fig. 9. Spin-down age vs period for the pulsars whose profiles are not 
affected by scattering in HBAs. Here the shape- and colour-coding rep¬ 
resent the difference of the number of peaks for each pulsar between 
HBAs and L-band relative to the HBAs. In red (triangles down) are the 
pulsars for which the number of peaks becomes smaller at HBAs than 
it is in L-band, in blue (triangles up) are the pulsars for which the num¬ 
ber of peaks becomes higher at HBAs than it is in L-band, and in green 
(triangles right) are the pulsars for which the number of peaks does not 
change. At the top and on the right-hand sides the stacked histograms 
represent the trends of the number of peaks as a function of period and 
age of the pulsar, respectively. 


The DM represents the integrated column density of free 
electrons between the pulsar and the observer. It produces a time 
delay in the signal, between the observing frequency and infinite 
frequency, that can be approximated as 


AfflM - 


r DM 1 
cm“-^pc 


2.41 ■ 10 - 4[^]2 


s. 


( 12 ) 


This approximation is valid if the plasma is tenuous and thus 
collisionless and if the observing frequency is much greater than 
the plasma frequency and the electron gyrofrequency. Signals at 
different frequencies will be delayed, with the lowest frequencies 
being delayed the most. These delays can change on timesca les 
of some years, up to 10“^ cm“^ pc (see e.g. lKeith et al.ll^oT^ . 

When aligning the profiles absolutely as we did, connect¬ 
ing the reference point of the profile to the reference epoch 
of the ephemeris, the DM delays had to be taken into ac¬ 
count and all reference times converted to the correspond- 
ing times at infinite freq uency to correct for dispersive delay. 
iHankins & Ricket^ (119861) described a way to measure DM vari¬ 
ations bas ed on the alignment o f pul sar profiles at differ e nt fre- 
quencies. IHankins et al.l (Il99ll) and IHankins & RankinI (1201 Ol) 
followed this method using increasingly higher resolution pro¬ 
files. They absolutely aligned the profiles spanning, where pos¬ 
sible, all the octaves of radio frequency. Leaving the DM as the 
only free parameter, they identified a reference point in the pro¬ 
file and adjusted the DM value to compensate for the remaining 
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misalignment of the profiles. The value of the best DM was ob¬ 
tained this way, its accuracy strongly depending on the lower 
frequency that can be used and on the precision with which the 
time difference of the misalignment can be measured. 

All the pulsars from our sample were aligned by refold¬ 
ing their profiles at all frequencies using the same ephemeris. 
The DM that was adopted for the alignment was the one ob¬ 
tained from LOFAR HBA data: a first folding was performed 
with prepfold using the new ephemeris created from the Lovell 
data (see Sect. O, but allowing for a search over DM values, 
and then the best DM obtained from this search was included 
in the ephemeris and the profile was dedeispersed once more, 
without any search option. The same was done for the LBA and 
the P- and L-band observations. Figure B1 shows that the multi¬ 
frequency profiles are aligned in most cases. There were some 
cases, discussed in Sect. |3l where the alignment is visibly in¬ 
correct or where we had to apply an extra adjustment to DM to 
compensate for a visible offset between the profiles (see Table 

D. 

One possible cause for the misalignment is the evolution of 
the profile across the frequencies, so that it is not possible to 
easily identify a fiducial point in the profile. In addition, pulsars 
strongly affected by scattering will not only have an exponential 
scattering tail, but also an absolute delay. Our profile alignment 
seemed to also be affected by the significantly different observ¬ 
ing epochs (e.g. the WSRT observations were performed more 
than ten years before LOFAR observations). On one hand, this 
meant that we had to adopt a timing solution spanning a long 
period of time, where timing noise and other effects can become 
substantial. On the other hand, with observations so far apart, we 
might also be probing gradients in electron density. 

In our case it is not yet possible to perform a system atic study 
of these variations, as was done bv iKeith et al.l (l2013h . Nonethe¬ 
less, LOFAR data can provide a new wealth of DM measure¬ 
ments to be compared with previous observations to map the 
evolution of the interstellar dispersion with time. A first con¬ 
clusion that can be drawn from this, simply by comparing the 
DM values in Table lB.fi is that there is no significant indication 
of DMs being systematically higher at low freque ncies (at least 
at our measurement precision). Some authors (e.g. lShit^ll983l) 
have postulated “superdispersion” due to the sweepback of field 
lines in the pulsar magnetosphere, which would be responsible 
for lower dispersion delays at low frequencies and could create 
an observed profile misalignment over a wide observing band. 
We find that the ratio DMeuA/DMeph ranges from 0.97 to 1.06 
(see Table EUi, thus differing by a significant percentage in 
some cases, but in both directions, thus not favouring superdis- 
persion. This agree s with previous findings usin g LOFAR data 
( [H^all et al.ll201^ and previous measurements (IHankins et ^ 

[IS 

5.3. Some examples of unexpected profile evolution 

It is not in the scope of this initial paper to enter into much detail 
about the profile evolution of specific pulsars. These will be the 
subject of future dedicated work. There are, however, several in¬ 
teresting cases of pulse evolution that are worth pointing out at 
this time. 

Figures [3] and [T] showed that there are some cases (e.g. 
B1541-t-09, B1821-H05, and B1822-09, B2224-H65) where the 
width of the profile is observed to increase with increasing ob¬ 
serving frequency. If we compare these results with the single 
profiles in Fig. lO we notice that in these cases this is caused 
by new peaks appearing in the profile at higher frequencies (e.g. 


the well-known ‘precursor’ in the case of B1822-09). This is 
not common in the standard core and cones geometries, even 
though it is predicted that new components can come into view if 
our viewing angle changes with increasing observing frequency, 
thus allowing us to see deeper into the beam. While this would 
explain some of the cases, in some others the profile evolution 
can hardly be ascribed to a symmetric core and conal structure 
(see e.g. B0355 -h54, B0450-h 55, B 1831-04, and B1857-26). 

These narrowing profiles at low frequen cies might be in¬ 
terpreted as evidenc e for fan beam models (iMic hell 198* ^. In 
fan beam models by iDvks et ^ (l2010h : iDvks & RudakI (1201 2l 
12015h . for instance, the emission comes from elongated broad¬ 
band streams that follow the magnetic field lines. The model, 
based on the cut angle at which the line of sight crosses the beam, 
can explain the lack of RFM, for example, if the stream is very 
narrow (also the case for millisecond pulsars), and it can also ex¬ 
plain the ‘inverse’ RFM if there is spectral non-uniformity along 
the azimuthal direction of the beam through which our line of 
sight c uts. The fan-beam formulation proposed by I Wang et ^ 
(l2014h can even explain ‘regular’ RFM by assuming that a fan 
beam composed of a (small) number of sub-beams will produce 
a so-called ‘limb-darkening pattern’ that is caused by the de¬ 
crease in intensity with beam radius of the emission at higher 
altitudes because the emission is farther from the magnetic pole. 
Their model, based on observations and simulations, predicts 
that the non-circularly bound beam (different in this from the 
beam predicted by the narrow-band models) can depart from the 
relation w oc F2 (where w is the measured width of the profile, 
see Sect. 14. 111 . 

IChen & Wangl(l2014h . who recently analysed the pulse width 
evolution with frequency of 150 pulsars from the EPN database, 
reached a similar conclusion: the emission must be broad-band, 
and the observed behaviour of width at different frequencies 
is caused by spectral changes along the flux tube. They found 
that if the spectral index variation along pulse phase is symmet¬ 
ric, there can be either canonical RFM or anti-RFM, while in 
cases where there are substantial deviations from the symmet¬ 
ric case, then there can be the non-monotonic trends of wio with 
V, which we also observed. This is supported theoretically by 
the p article-in-cell sim ulations of pair production in the vacuum 
gaps (lTimokhinll20 1 Oh . which predict that the secondary plasma 
does not necessarily have a monotonic momentum spectrum. 

These results, and in particular the fact that a wide stream is 
expected to produce spectral variations longitudinally in the pro¬ 
files, could also explain the observed behaviour of the peak ratios 
shown in Fig. |6] Alternatively, the peak ratio changing with ob¬ 
serving frequency and, in particular, its changing sign, might be 
related to the frequency dependance of the two mod es of polar¬ 
isatio n (ordinary, ‘O’, and extraordinary, ‘X’) (e.g.. lSmits et al.l 
l2006h and to that they might be differently dominant in different 
peaks. While it is not possible to give a comprehensive analysis 
of the phenomenon here, our s tudies on the polaris ed emission 
from pulsars with LOFAR (see lNoutsos et al.ll2015h address the 
questions related to the orthogonally polarised modes and the 
related jumps in the polarisation angle. 

While in this section we have discussed some unexpected 
profile evolution, it remains the case that most pulsar profiles 
(cf. Sects. l4~n and l5.LTT i are well described by RFM. 


6. Summary and future work 

We have presented the profiles of the first 100 pulsars observed 
by LOFAR in the frequency range 119-167 MHz. Twenty-six 
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of them were also detected with 57 min integrations using LO¬ 
FAR in the interval 15-63 MHz. All the LOFAR profiles pre¬ 
sented in this work will be made available through the EPN 
databases LOFAR observations were compared with archival 
WSRT or Lovell observations in P- and L-band, after first fold¬ 
ing and aligning all profiles using an ephemeris spanning the full 
range of the observations. The rotational and derived parameters 
are presented in Table IB. Il Two values of DM were presented as 
well: one obtained by the best timing fit and one from the best fit 
of LOFAR data. For each pulsar we aligned the profiles at dif¬ 
ferent frequencies in absolute phase, using the latter DM value. 
The 100 profiles are presented in Fig. lO 

Each profile of every pulsar was described using a m ulti- 
Gaussian fit following the approach of iKramer et al.l 119941) . so 
that in general more components were needed to fit the pro- 
files than evide nt at first glance or traditionally considered (e.g. 
lRankiniri983bll . The results of the Gaussian fit (the measure¬ 
ment of the widths at half and at 10% of the maximum of the 
profile and the spectral index of their evolution with observing 
frequency) are reported in Table IB. 21 Using the components’ 
widths, we calculated the ratio of the peaks for pulsars with dou¬ 
ble or multiple peaks (the two most prominent, in the case of 
multiple peaks pulsars). We concluded that the ratio of the main 
peak to the second peak does not follow a unique trend, although 
we note that in most cases the dominant peak alternates with 
changing observing frequency. Using ivw , we followed the evo¬ 
lution of the width of the full profile with observing frequency. 
We concluded that while our average spectral index is compati¬ 
ble with no evolution of the pulse width, the distribution of the 
values is quite large and compatible with the presence of differ¬ 
ent behaviours for different pulsars, for example, based on inne r 
or outer cone emission, as discussed bv iMitra & RankinI (l2002l) . 
or on different propagation modes in the m agnetosphere (e.g. 
iBeskin et al.lll988tlBeskin & PhiliDDOvll2012h . 

Euture work will be needed, and is in progress, to add more 
elements to complete this puzzle. Parallel to this work, a similar 
one is being conducted on the evol ution of the profiles of mil¬ 
lisecond pulsars at low frequencies (iKondratiev et al.ir2015h . and 
the spectral behaviour of the slow and recycled pulsars has been 
analysed (Hassall et al. in prep.). Complementary to our work is 
the st udy of the polarisation properties of pulsars (iNoutsos et alJ 
l20L5h : the study of the polarisation properties can give con¬ 
straints on the geometry of the pulsar emission and therefore 
on its height and on the intrinsic opening angle of the beam 
of the emission. Additionally, polarisation can help distinguish 
between orthogonal polarisation modes and therefore determine 
whether the observed widening of the profiles is caused by bire¬ 
fringence. 

To better constrain the width of the pulse, it is also impor¬ 
tant to be able to deconvolve the scattering tail, which in some 
cases becomes dominant at low frequencies, from the intrinsic 
width of the profile. Studies on the characterisation of scattering 
at low frequen cies and modelling of the scattering tail are be¬ 
ing conducted (lArchibald et alJl20l4 Zagkouris et al. in prep.). 
At the same time, the effects of the ISM on LOEAR profiles are 
being used to create an TSM weather’ database, where the DM 
variations, independently measured with LOEAR (Verbiest et al. 
in prep.), can be used for hig h-precision timin g measurements 
from pulsar timing arrays (see lKeith et al.ll2()T^ . 

Einally, the observations presented in this work were from 
commissioning LOEAR data; as mentioned, LOEAR is currently 
performing the LOTAAS all-sky pulsar survey. When com- 

http://www.epta.eu.org/epndb/ 


pleted, it will provide an extraordinary database of 1.5 h cover¬ 
age of the whole northern sky with 0.49 ms time resolution and 
down to a flux ~ 6 mJy at 135 MHz. At present we are able 
to use the full LOEAR core instead of only the ‘Superterp’ and 
can cover a full 80 MHz bandwidth contiguously. Moreover, fur¬ 
ther improvements will include coherent dedispersion of the data 
and will enable single pulse studies to probe the ‘instantaneous’ 
magnetosphere. 
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Appendix A: Determination of puise widths 

As discussed in the Sects. [^and l4. II the widths of the full profile that are to be used in the subsequent calculations were determined 
from the fit of Gaussian components to the profile shape. However, a number of alternative methods were also tested and, in one 
case, used to cross check the widths obtained from the Gaussian fits. Here we give a brief description of each of them. 

We calculated the effective width, Wes, as the integrated pulsed flux divided by the peak flux; this is represented by the cyan- 
shaded areas in Fig. lA. 1 I centred at the main peak. This width metric typically does not represent all the profile components well, as 
is observed in the P- and L-band profiles. The method underestimates the width of the profile for multi-peaked profiles, which is 
the majority of profiles in this work, mainly characterising the width of the main peak in these cases. 

The total-power width Wpow is calculated as the phase interval that includes 90% of the total pulse energy in the cumulative flux 
distribution of the normalised profile (see lower plots in Fig. lA.lb . We selected the phase included in the interval between 5% of the 
total flux and 95% of the total flux (indicated by the intersection of the summed profile with the horizontal lines in panel (b) of Fig. 
lA. lb . This distribution should start at 0 and increase monotomically to 1 if no noise is present (see the noise-free distribution on the 
right side of Fig. lATb . When noise is present (left-hand side of Fig. lATI) . the sum does not grow smoothly, and sometimes negative 
terms in the off-pulse region add up quite substantially. In the real profiles (left-hand side of Fig. lATb . the horizontal lines should 
all overlap at 0.05 and 0.95, like they do in the case of the noise-free profiles (right-hand side of Fig. lA.lb . but they differ in some 
cases as they take into account that the highest and lowest value of a noisy cumulative distribution can be negative or greater than 
1. As can be seen from the green (LBA) and magenta (P-band) curves, when the profile is noisier, this measurement is less stable 
because the cumulative distribution oscillates more, and the width can be consistently overestimated (in these cases the left-hand 
dashed line precedes the phase range shown in the Fig.). This method is represented by the dotted line in the upper plot of Fig. IA.il 
which, in the case of the LBA and P-band profiles, starts at phase 0. 

The full-width at 10% of the maximum, Wop, represents the width of the full profile (including noise) at the 10% level of the 
outer components, including the full on-pulse region: it is coincident with the solid vertical lines in Fig. IA.il It was calculated as 
the on-pulse region with a flux above 10% of the main peak in the baseline-subtracted profile. 
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Spin Phase 
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Fig. A.l. PSR B1112+50: example of how the widths of the full profile were calculated using the four methods described in the text. On the 
left-hand side is the real profile. On the right-hand side, for comparison, the noise-free profile obtained from the Gaussian fit. In panel (a), the 
shaded cyan area represents Wcb- The dashed vertical lines represent Wpow, calculated as shown in panel (b): the cumulative distribution of the 
normalised profile. The solid line demarcates wio as obtained from the Gaussian smoothed profiles (which are shown in panel (a), on the right). 
The measurement of ww is consistent and overlapping with iPop (see text for details). 
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Appendix B: Multi-frequency profiles and tables of the derived profile properties. 


Fig. B.l. Pulse profiles of the 100 pulsars observed by LOFAR, showing from bottom to top LB A (green) and HBA (blue) profiles from LOFAR 
observations, 350 MFIz (magenta) and 1400 MFIz (red) profiles from WSRT and Lovell observations. We caution that the moding pulsars, for 
which only one mode of pulsation is represented by our observation, might not be the same for the profiles at different frequencies. The profile 
phase is zoomed-in on the interval 0.25 - 0.75 in all cases except for the pulsars with interpulses, where the full pulse phase is shown, and the 
profile is rotated by 0.25 in phase to shift the interpulse from phase 0. The star next to the band name indicates that the alignment for that band 
was made manually (see text for details). 
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Fig IB.ll continued 
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Fig IB.ll continued 
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Fig. B.4. * 


Fig IB.ll continued 
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Fig IB.ll continued 
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Fig. B.6. * 
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Fig. B.7. Zoom-in on the 0.0 - 0.5 phase longitude of the profile for the five pulsars in our sample thaf have inferpulses. The plots are scaled so 
that each profile has the same r.m.s. and are renormalised to the peak of the interpulse at each frequency. In all cases LBA profiles were removed, 
as their signal-to-noise ratio and their channel smearing did not enable any estimate of the interpulse. The reference phase of the main pulse has 
been shifted by 0.75, relative to the main plots, to place the interpulse at the centre of the plot. 



































Table B.l. List of the 100 pulsars and their rotational and derived parameters. 


PSR Name 

P 

[s] 

P 

[s s-‘] 

Epoch 

[MJD] 

LBA Epoch 
[MJD] 

HBA Epoch 
[MJD] 

[cm“'^pc] 

77Mhba 

[cm^^pc] 

logT 

[yr] 

logB 

[G] 

logE 

[s-'erg] 

Notes 

BOO11+47 

1.241 

5.647e-16 

55316 


55858 

30.67(7) ['i 

30.3(5) 

7.54 

11.93 

31.07 

WES"' 

B0031-07 

0.943 

4.081e-16 

55857 


56630 

11.28(8) ['1 

10.9(4) 

7.56 

11.80 

31.28 

WES"' 

J0051+0423 

0.355 

7.047e-18 

49800 


56318 

13.9(1) P] 

13.93(1) 

8.90 

10.70 

30.79 

LOTAAS* 

B0059+65 

1.679 

5.970e-15 

56475 


55859 

65.85(2) ['1 

65.8(8) 

6.65 

12.51 

31.70 

EPN"' 

BO105+65 

1.284 

1.292e-14 

54578 


55859 

30.6(5) 

30.6(2) 

6.20 

12.62 

32.38 

P- 2004" 

BOl14+58 

0.101 

5.846e-15 

55829 


55859 

49.462(4) 

49.42(1) 

5.44 

11.89 

35.34 

EPN"! 

B0136+57 

0.272 

1.069e-14 

54589 


55859 

73.776(6) 

73.83(1) 

5.61 

12.24 

34.32 

P-, L- 2004* 

B0138+59 

1.223 

3.910e-16 

54266 

55924 

55859 

34.74(1) ['1 

34.92(2) 

7.70 

11.84 

30.93 

P- 2004" 

B0301+19 

1.388 

1.295e-15 

55743 


55859 

15.851(9) 

15.65(4) 

7.23 

12.13 

31.28 

WES"' 

B0320+39 

3.032 

6.358e-16 

54125 

55924 

55859 

25.99(3) 

26.19(1) 

7.88 

12.15 

29.95 

P- 2002" 

B0329+54 

0.715 

2.046e-15 

53850 

55859 

55859 

26.7941(1) [31 

26.7653(< 1) 

6.74 

12.09 

32.35 

P- 2000" 

B0331+45 

0.269 

7.331e-18 

55899 


55859 

47.149(3) I'l 

47.15(1) 

8.76 

10.65 

31.17 

EPN"' 

B0339+53 

1.934 

1.342e-14 

55804 


55859 

67.28(6) ['1 

67.5(2) 

6.36 

12.71 

31.86 

WES"' 

B0355+54 

0.156 

4.395e-15 

53813 


55859 

57.1226(3) W 

57.14(1) 

5.75 

11.92 

34.66 

P- 2003" 

B0402+61 

0.595 

5.576e-15 

56097 


55859 

65.303(7) [11 

65.39(3) 

6.23 

12.27 

33.02 


B0410+69 

0.391 

7.656e-17 

54418 


55859 

27.471(4) [11 

27.44(1) 

7.91 

11.24 

31.70 

P- 2004" 

B0447-12 

0.438 

1.028e-16 

55203 


56568 

37.04(1) ['1 

37.01(1) 

7.83 

11.33 

31.68 

P- 2003" 

B0450+55 

0.341 

2.377e-15 

54139 

55924 

55861 

14.483(7) ['1 

14.59(7) 

6.36 

11.96 

33.38 

P-, L- 2004" 

B0450-18 

0.549 

5.755e-15 

55203 


56603 

39.903(3)1'! 

39.91(2) 

6.18 

12.25 

33.14 

P- 2004" 

B0523+11 

0.354 

7.362e-17 

55202 


56568 

79.423(3) ['! 

79.43(4) 

7.88 

11.21 

31.81 

P- 2006" 

B0525+21 

3.746 

4.004e-14 

54457 

55924 

56596 

50.96(2) ['1 

50.87(2) 

6.17 

13.09 

31.48 

P- 2003, L- 2004" 

B0540+23 

0.246 

1.542e-14 

54748 


55859 

77.631(2)1'! 

77.7(9) 

5.40 

12.29 

34.61 

P- 2006" 

B0609+37 

0.298 

5.980e-17 

54590 


55859 

27.136(4)1'! 

27.15(3) 

7.90 

11.13 

31.95 

P- 2006" 

B0611+22 

0.335 

5.911e-14 

55604 


56681 

96.94(7) ['! 

96.92(1) 

4.95 

12.65 

34.79 

WES"' 

B0626+24 

0.477 

1.985e-15 

55071 


55859 

84.19(4)1'! 

84.18(3) 

6.58 

11.99 

32.86 

P- 2004" 

B0643+80 

1.214 

3.799e-15 

54855 


55859 

33.33(2) ['! 

33.3(1) 

6.70 

12.34 

31.92 

P- 2004" 

B0809+74 

1.292 

1.660e-16 

53508 

56046 

55916 

5.733(1)13! 

5.75(3) 

8.09 

11.67 

30.48 

WES"' 

B0818-13 

1.238 

2.105e-15 

55507 


56603 

40.949(3) 1'! 

40.98(1) 

6.97 

12.21 

31.64 

WES"' 

B0820+02 

0.865 

1.047e-16 

55508 


56301 

23.717(6)1'! 

23.76(7) 

8.12 

11.48 

30.81 

WES"' 

B0823+26 

0.531 

1.679e-15 

55419 

55889 

55966 

19.464(4) 1'! 

19.4711(< 1) 

6.70 

11.98 

32.65 

WES"' 

B0834+06 

1.274 

6.799e-15 

55507 

55883 

55881 

12.891(6)1'! 

12.8735(< 1) 

6.47 

12.47 

32.11 

WES"' 

B0906-17 

0.402 

6.693e-16 

54887 


56603 

15.885(3)1'! 

15.87(1) 

6.98 

11.72 

32.61 

WES"' 

B0917+63 

1.568 

3.608e-15 

53883 

55925 

55860 

13.21(2)1'! 

13.14(2) 

6.84 

12.38 

31.57 

EPN"! 

B0919+06 

0.431 

1.372e-14 

55870 


56315 

27.256(6) 1'! 

27.3(1) 

5.70 

12.39 

33.83 

WES"! 

B0943+10 

1.098 

3.546e-15 

54226 

55925 

55869 

15.2(5)1'! 

15.32(1) 

6.69 

12.30 

32.02 

EPN"! 

B0950+08 

0.253 

2.165e-16 

51541 

55889 

56603 

2.961(3)1'! 

2.9696(< 1) 

7.27 

11.37 

32.72 

EPN"! 

B1112+50 

1.656 

2.493e-15 

54121 

55925 

56233 

9.202(8) 1'! 

9.18(1) 

7.02 

12.31 

31.34 

WES"! 

B1133+16 

1.188 

3.733e-15 

56235 

55925 

55916 

4.8611(1)13! 

4.8262(< 1) 

6.70 

12.33 

31.94 

WES"' 

B1237+25 

1.382 

9.596e-16 

46498 

55925 

55860 

9.254(6) I'l 

9.2538(< 1) 

7.36 

12.07 

31.16 

WES"' 

J1238+21 

1.119 

1.445e-15 

55507 


55860 

17(2) ['! 

17.97(3) 

7.09 

12.11 

31.61 
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Table B.l. Continued. 


PSR Name 

P 

[s] 

P 

[s s^‘] 

Epoch 

[MJD] 

LBA Epoch 
[MJD] 

HBA Epoch 
[MJD] 

£>Meph“ 

[cm^'^pc] 

77Mhba 

[cm“5pc] 

logT 

[yr] 

logB 

[G] 

logE 

[s^'erg] 

Notes 

J1313+0931 

0.849 

8.000e-16 

50984 


55860 

12 .0(1) 

12.05(2) 

7.23 

11.92 

31.71 


B1322+83 

0.670 

5.675e-16 

53763 


55860 

13.3(2) ['[ 

13.31(3) 

7.27 

11.80 

31.87 

WES'' 

B1508+55 

0.740 

5.007e-15 

55168 

55906 

55860 

19.59(2) ['[ 

19.6179(< 1) 

6.37 

12.29 

32.69 

WES'' 

B1530+27 

1.125 

7.791e-16 

56508 


55859 

14.76(2) ['[ 

14.7(2) 

7.36 

11.98 

31.33 

WES'' 

B1540-06 

0.709 

8.810e-16 

55962 

55892 

56605 

18.403(4) ['[ 

18.3788(< 1) 

7.11 

11.90 

31.99 

WES'' 

B1541+09 

0.748 

4.328e-16 

55508 

55892 

55859 

35.24(3) ['[ 

34.97(3) 

7.44 

11.76 

31.61 

WES'' 

B1604-00 

0.422 

3.063e-16 

54696 

55892 

56346 

10.687(5) ['[ 

10.69(1) 

7.34 

11.56 

32.21 

WES'' 

B1612+07 

1.207 

2.360e-15 

55508 


56354 

21.39(3) ['[ 

21.41(1) 

6.91 

12.23 

31.72 

WES'' 

J1627+1419 

0.491 

3.930e-16 

48957 


55859 

33.8(6) [5[ 

32.18(9) 

7.30 

11.65 

32.12 


B1633+24 

0.491 

1.192e-16 

55508 

56517 

55859 

24.33(4) ['[ 

24.27(1) 

7.81 

11.39 

31.60 

EPN'' 

B1642-03 

0.388 

1.770e-15 

55643 


56605 

35.738(3) [[[ 

35.7556(< 1) 

6.54 

11.92 

33.08 

WES'' 

J1645+1012 

0.411 

8.203e-17 

52814 


56276 

36.4(6) [51 

36.18(3) 

7.90 

11.27 

31.67 


J1652+2651 

0.916 

6.542e-16 

55507 


55860 

41.0(2) [2] 

40.79(2) 

7.35 

11.89 

31.53 

WES'' 

B1737+13 

0.803 

1.441e-15 

55507 


55870 

48.672(4) ['[ 

48.67(1) 

6.95 

12.04 

32.04 

WES'' 

B1749-28 

0.563 

8.000e-15 

54111 


56606 

50.39(8) ['[ 

50.38(1) 

6.05 

12.33 

33.25 

WES'' 

B1818-04 

0.598 

6.347e-15 

56445 


56606 

84.42(2) ['[ 

84.34(3) 

6.17 

12.29 

33.07 

WES'' 

B1821+05 

0.753 

2.267e-16 

55508 


56286 

66.777(3) [[[ 

66.79(2) 

7.72 

11.62 

31.32 

LOTAAS* 

B1822-09 

0.769 

5.249e-14 

56183 

55892 

56605 

19.36(4) ['[ 

19.36(2) 

5.37 

12.81 

33.66 

WES'' 

B1831-04 

0.290 

7.250e-17 

55508 


56605 

79.302(8) [[[ 

79.39(4) 

7.80 

11.17 

32.07 

WES'' 

B1839+09 

0.381 

1.089e-15 

55541 


55870 

49.145(8) ['[ 

49.1501(< 1) 

6.74 

11.81 

32.89 

WES'' 

B1839+56 

1.653 

1.495e-15 

55508 

55892 

55870 

26.71(1) ['1 

26.7540(< 1) 

7.24 

12.20 

31.12 

WES'' 

B1842+14 

0.375 

1.873e-15 

55817 

55864 

55870 

41.512(4) ['[ 

41.48(1) 

6.50 

11.93 

33.15 

WES'' 

B1848+12 

1.205 

1.137e-14 

55835 


56276 

70.61(2) ['[ 

70.6(8) 

6.23 

12.57 

32.41 

LOTAAS'' ,WES'' 

B1848+13 

0.346 

1.493e-15 

55508 


56276 

60.149(8) ['[ 

60.1(1) 

6.56 

11.86 

33.15 

LOTAAS'’ ,EPN'' 

B1857-26 

0.612 

2.050e-16 

55508 


55870 

37.993(5) [[[ 

37.92(3) 

7.68 

11.55 

31.55 

WES'' 

B1905+39 

1.236 

5.408e-16 

55507 


55870 

30.96(3) ['[ 

30.98(5) 

7.56 

11.92 

31.05 

WES'' 

B1907+00 

1.017 

5.513e-15 

55770 


56354 

112.776(6) [[[ 

113(1) 

6.47 

12.38 

32.32 

LOTAAS'’ ,WES'' 

B1907+02 

0.990 

5.492e-15 

53773 


56339 

171.722(9) [[[ 

172(1) 

6.46 

12.37 

32.35 

LOTAAS* ,WES'' 

B1907+10 

0.284 

2.644e-15 

55729 


55870 

149.756(4) [[[ 

149.9(1) 

6.23 

11.94 

33.66 

WES'' 

B1911-04 

0.826 

4.069e-15 

55705 


56354 

89.38(1) ['[ 

89.38(3) 

6.51 

12.27 

32.46 

LOTAAS'' ,WES'' 

B1914+09 

0.270 

2.518e-15 

55931 


55870 

60.981(6) ['[ 

61.01(6) 

6.23 

11.92 

33.70 

WES'' 

B1915+13 

0.195 

7.194e-15 

55745 


55870 

94.54(4) ['[ 

94.65(3) 

5.63 

12.08 

34.59 

WES'' 

B1917+00 

1.272 

7.675e-15 

55695 


56354 

90.27(2) ['[ 

90(1) 

6.42 

12.50 

32.17 

LOTAAS'’ ,WES'' 

B1918+26 

0.786 

3.414e-17 

54142 


55870 

27.68(9) ['[ 

27.7(1) 

8.56 

11.22 

30.44 

EPN'' 

B1919+21 

1.337 

1.349e-15 

55969 

55860 

55859 

12.4632(1) [5[ 

12.4453(< 1) 

7.20 

12.13 

31.35 

WES'' 

B1920+21 

1.078 

8.169e-15 

55507 


55859 

217.08(2) [[[ 

217.0(1) 

6.32 

12.48 

32.41 

WES'' 

B1923+04 

1.074 

2.459e-15 

55445 


56339 

102.24(2) [[[ 

102.2(9) 

6.84 

12.22 

31.89 

LOTAAS'’ ,WES'' 

B1929+10 

0.227 

1.157e-15 

55830 

55913 

55859 

3.178(4) ['[ 

3.18(1) 

6.49 

11.71 

33.59 

WES'' 

B1953+50 

0.519 

1.374e-15 

55508 

55864 

55858 

31.974(3) ['[ 

31.98(1) 

6.78 

11.93 

32.59 

EPN'' 
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Table B.l. Continued. 


PSR Name 

P 

[s] 

P 

[s s^‘] 

Epoch 

[MJD] 

LBA Epoch 
[MJD] 

HBA Epoch 
[MJD] 

£>Meph“ 

[cm^'^pc] 

77Mhba 

[cm^^pc] 

logT 

[yr] 

logB 

[G] 

logE 

[Y'erg] 

Notes 

B2021H-51 

0.529 

3.066e-15 

55889 


55858 

22.634(6) '‘1 

22.55(1) 

6.44 

12.11 

32.91 

WES"' 

B2022H-50 

0.373 

2.515e-15 

55507 


55859 

33.036(3) 

33.0(2) 

6.37 

11.99 

33.28 

WES"' 

B2043-04 

1.547 

1.471e-15 

55311 


56324 

35.81(1) ['1 

36(1) 

7.22 

12.18 

31.20 

WES"' 

B2044H-15 

1.138 

1.823e-16 

55507 


56293 

39.86(1) 

39.8(1) 

8.00 

11.66 

30.69 

LOTAAS* ,P-2004" 

B2106H-44 

0.415 

8.606e-17 

55855 


56602 

139.83(1) 

140(1) 

7.88 

11.28 

31.68 

WES"' 

B2110H-27 

1.203 

2.623e-15 

55424 


56602 

25.119(4) 

25.1(1) 

6.86 

12.25 

31.77 

WES"' 

B2113H-14 

0.440 

2.895e-16 

55696 


55859 

56.15(7) 

56.2(1) 

7.38 

11.56 

32.13 

WES"' 

B2122H-13 

0.694 

7.677e-16 

55507 


56332 

30.1(1) 

30.3(2) 

7.16 

11.87 

31.96 

LOTAAS* ,EPN"' 

B2148H-63 

0.380 

1.707e-16 

55769 


55858 

129.704(1) 

129.7(2) 

7.55 

11.41 

32.09 

WES"' 

B2154H-40 

1.525 

3.439e-15 

55334 


55858 

70.86(1) 

71.14(6) 

6.85 

12.37 

31.58 

WES"' 

B2217H-47 

0.538 

2.767e-15 

56228 

55860 

55858 

43.49(1) 

43.4811(< 1) 

6.49 

12.09 

32.84 

WES"' 

B2224H-65 

0.683 

9.647e-15 

55964 


55858 

36.08(9) 

36.33(5) 

6.05 

12.41 

33.08 

WES"' 

B2227H-61 

0.443 

2.255e-15 

55859 


55858 

124.62(2) 

124.7(1) 

6.49 

12.00 

33.01 

WES"' 

J2248-0101 

0.477 

6.596e-16 

56405 


56304 

29.06(3) 

29.1(4) 

7.06 

11.75 

32.38 

LOTAAS* ,WES^ 

J2253-H1516 

0.792 

6.653e-17 

49079 


56354 

29.182(9) ® 

29(1) 

8.28 

11.37 

30.72 

LOTAAS* 

B2255H-58 

0.368 

5.749e-15 

55671 


55858 

151.082(6) 

151.2(1) 

6.01 

12.17 

33.66 

WES"' 

B2303H-30 

1.576 

2.893e-15 

54584 

55892 

55858 

49.56(2) ['1 

49.58(5) 

6.94 

12.33 

31.47 

WES"' 

B2306H-55 

0.475 

1.995e-16 

55384 


55858 

46.535(8) 

46.54(2) 

7.58 

11.49 

31.87 

WES"' 

B2310H-42 

0.349 

1.123e-16 

54363 

55860 

56606 

17.274(1) 

17.27(1) 

7.69 

11.30 

32.02 

EPN"' 

B2315H-21 

1.445 

1.047e-15 

56407 


55858 

20.874(7) 

20.87(6) 

7.34 

12.10 

31.14 

P- 2004" 

B2334H-61 

0.495 

1.933e-13 

55508 


55858 

58.41(2) 

58.43(3) 

4.61 

13.00 

34.80 

WES"' 


Notes. For each DM measurement and its error we quote the reference paper: 1, 2, 3, 4, 5, 6, 7, and 8. Observations from the LOTAAS survey dCoenen et al.ll201^ . Year of the P- 
and/or L-hand data from WSRT observations. L-band data were obtained from Lovell observations in all but the three cases where the year of t he WSRT observations i s indicated. All the L-band 
observations using the Lovell telescope were performed between October 2011 and November 2013. Where P-band data were available from lWeltevrede et al.l (l2005h . but no absolute reference 
was available to align the data with LOFAR data, reference is made to their paper as WES and the corresponding profile in Fie. lB.ll was aligned by eye. For completeness, where the pulsar was 
not in WES, then EPN indicates that the profile in P-band was obtained from the EPN pulsar database and was also only aligned by eye. In this latter case, the observations are not later than 2002. 
In the cases where nothing is specified, for the P-band the profile is either missing in all the databases used here or the Lovell data were used and the profile aligned with the standard procedure as 
described in the text. 

Refer ence s. fllFIobbs et al. 

(Il987h : rSi lCamilo&Nicd 


j200^; f2i lSaver et al.l ( Il997l) : Gl lHassall et al.l ( l2012h : 14l lJanssen & Stappersl ( 120061) : fSl lLommen et al.l ( l2000h : (hl lLewandowski et al.l ( |2004|) : (7l lFtamilton & L^^ 
199.5h . 
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Table B.2. Width of the full profile and duty cycle. For each pulsar (Col. 1) at each observing band (Col. 2), we list in Col. 3 the full width at 
half maximum, wso, in Col. 4 the width of the full profile calculated as wio (see text for details and Fig. lA. lb . an d in Col. 5 is the duty cycle of the 
pulse, wiolP, indicated as percent age of the total profile. In the last two columns we list the spectral indices 6 of the evolution of w^o and wiq with 
observing frequency, modelled as w(v) oc i/. In all cases the statistical error is quoted in parentheses on the last digit. Pulsars with an interpulse 
are marked with the notation (IP) next to their names. Where scattering affects the HB A measurement, no spectrum was produced (see Fig.[3]for 
reference). 


PSR Name 

Band 

W50 

ww 

WiolP 

^50 

<5io 



[deg] 

[deg] 

% 




BOO11+47 

BOOS1-07 

J0051+0423 
B0059+65 

BO105+65 

BOl14+58 

B0136+57 

B0138+59 

B0301+19 

B0320+39 

B0329+54 

BOSS1+45 

B0339+53 


HBA 

22(1) 

43(3) 

P 

21(4) 

38(2) 

L 

15.8(3) 

44.7(6) 

HBA 

26(2) 

45(3) 

P 

21.1(4) 

39.9(4) 

L 

11.3(3) 

34.5(3) 

HBA 

15.5(3) 

36.6(5) 

HBA 

2(1) 

26(1) 

P 

18.3(8) 

24(2) 

L 

16.2(3) 

19.4(3) 

HBA 

8.4(6) 

16(1) 

P 

4(1) 

15(1) 

L 

7.0(3) 

16.2(4) 

HBA 

scattered 


P 

10(2) 

31(2) 

L 

10(1) 

16.5(7) 

HBA 

scattered 


P 

4(1) 

9(1) 

L 

4(1) 

9(1) 

LBA 

5.6(7) 

50(4) 

HBA 

5.6(3) 

32(1) 

P 

8(1) 

31(1) 

L 

12.3(3) 

32.7(5) 

HBA 

17.6(6) 

22.5(6) 

P 

15.3(6) 

15.5(6) 

L 

10.9(3) 

14.8(3) 

LBA 

8(2) 

15(2) 

HBA 

4.6(3) 

9.8(3) 

P 

4(1) 

8(1) 

L 

6.7(8) 

11.6(3) 

LBA 

scattered 


HBA 

8.8(3) 

32.7(3) 

P 

2(1) 

29(1) 

L 

1.4(3) 

26.8(3) 

HBA 

5(2) 

19(2) 

P 

7.6(9) 

17(1) 

L 

1.8(3) 

9.9(4) 

HBA 

7.7(6) 

16.2(7) 

P 

7(1) 

11(1) 

L 

5.3(5) 

8.4(7) 


-0.16(8) 0.04(7) 


12.1(9) 

10.6(7) 

12.4(1) 

-0.46(6) 

-0.12(3) 

12.5(8) 

11.1(1) 

9.6(1) 

-(-) 

-(-) 

10.2(1) 

-0.08(9) 

-0.15(7) 

7.4(5) 

6.8(6) 

5.38(9) 

-0.07(9) 

0.01(9) 

4.5(4) 

4.3(3) 

4.5(1) 

-(-) 

-(-) 

8.8(6) 

4.6(2) 

-(-) 

-(-) 

2.7(3) 

2.7(3) 

0.3(6) 

-0.04(3) 

14(1) 

9.1(3) 

8.6(3) 

9.1(1) 

-0.21(5) 

-0.18(4) 

6.3(1) 

4.3(1) 

4.11(9) 

0.1(1) 

0.05(4) 

4.2(6) 

2.73(9) 

2.4(3) 

3.23(9) 

0.1(1) 

-0.09(2) 

9.08(9) 

8.2(3) 

7.45(9) 

-0.9(3) 

-0.4(1) 

5.5(7) 

4.7(3) 

2.7(1) 

-0.2(1) 

-0.3(1) 

4.5(2) 

3.1(3) 

2.3(2) 

0.1(4) 

0.09(3) 


B0355+54 
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Table B.2. Continued. 


PSR Name 

Band 

W50 

ww 

wwlP 

^50 

<5io 



[deg] 

[deg] 

% 




HBA 

13.0(4) 

30.2(8) 

8.4(2) 




P 

5(1) 

18(1) 

5.1(3) 




L 

16.2(3) 

35.5(3) 

9.87(9) 



B0402+61 

HBA 

16(1) 

28(1) 

7.8(3) 

-0.8(1) 

-0.23(5) 


P 

3.5(9) 

15.5(6) 

4.3(1) 




L 

2.8(3) 

14.8(3) 

4.11(9) 



B0410+69 

HBA 

5(2) 

19(2) 

5.5(7) 

-0.4(4) 

-0.4(1) 


P 

2(1) 

12(1) 

3.5(3) 




L 

2.1(3) 

8.4(8) 

2.3(2) 



B0447-12 

HBA 

2.8(2) 

25(2) 

7.0(7) 

0.8(1) 

0.1(1) 


P 

12(1) 

23(1) 

6.7(3) 




L 

20(1) 

27(1) 

7.7(5) 



B0450+55 

LBA 

5(2) 

39(4) 

10(1) 

0.1(1) 

-0.07(5) 


HBA 

25(4) 

39(2) 

10.9(7) 




P 

8(1) 

28(2) 

7.8(6) 




L 

17.9(3) 

31.0(3) 

8.6(9) 



B0450-18 

HBA 

21(1) 

32(2) 

9.0(7) 

-0.1(1) 

-0.1(1) 


P 

16(1) 

22(1) 

6.3(3) 




L 

16(1) 

22(1) 

6.3(3) 



B0523+11 

HBA 

scattered 



-(-) 

-(-) 


P 

13(1) 

21(1) 

5.9(3) 




L 

13.4(3) 

17.2(3) 

4.79(9) 



B0525+21 

LBA 

25(1) 

28(2) 

8.0(6) 

-0.17(6) 

-0.12(6) 


HBA 

21.3(1) 

25.3(1) 

7.03(4) 




P 

16(1) 

21(1) 

5.9(3) 




L 

14(1) 

19(1) 

5.5(3) 



B0540+23 

HBA 

scattered 



-(-) 

-(-) 


P 

7(1) 

23(1) 

6.7(3) 




L 

7.4(3) 

24.6(3) 

6.84(9) 



B0609+37 

HBA 

2(2) 

30(2) 

8.6(7) 

0.1(2) 

-0.23(8) 


P 

4.5(3) 

24(1) 

6.7(3) 




L 

4.9(3) 

17.9(4) 

5.0(1) 



B0611+22 

HBA 

scattered 



-(-) 

-(-) 


P 

5(1) 

14(1) 

3.9(4) 




L 

6.7(3) 

14.1(3) 

3.91(9) 



B0626+24 

HBA 

7(1) 

23(1) 

6.6(4) 

-(-) 

-(-) 


P 

5(1) 

14(1) 

3.9(3) 




L 

9.9(3) 

17.2(3) 

4.79(9) 



B0643+80 

HBA 

4.9(6) 

9(1) 

2.7(3) 

-0.3(1) 

-0.2(1) 


P 

1(1) 

11(1) 

3.1(3) 




L 

2.5(3) 

6.7(3) 

1.86(9) 



B0809+74 

LBA 

33.4(3) 

52.7(9) 

14.6(2) 

-0.32(1) 

-0.17(1) 


HBA 

12.3(3) 

29.5(3) 

8.2(1) 




P 

12.1(2) 

25.6(2) 

7.1(6) 




L 

14.4(3) 

26.4(3) 

7.33(9) 



B0818-13 

HBA 

4.9(6) 

10.5(6) 

2.9(1) 

0.0(1) 

-0.06(6) 
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PSR Name 

Band 

W50 

ww 

WlolP 

^50 

<5io 



[deg] 

[deg] 

% 




B0820+02 

B0823+26 (IP) 

B0834+06 

B0906-17 

B0917+63 

B0919+06 

B0943+10 

B0950+08 (IP) 

BI112+50 

B1133+16 

B1237+25 

J1238+21 

J1313+0931 


p 

6.2(1) 

9.8(1) 

L 

6.0(3) 

9.1(3) 

HBA 

12(1) 

19(1) 

P 

9.2(3) 

14.3(3) 

L 

7.7(3) 

13.0(3) 

LBA 

10.5(3) 

34(2) 

HBA 

4(1) 

12(1) 

P 

3.3(7) 

9.2(7) 

L 

2.8(3) 

6.0(3) 

LBA 

7.7(3) 

15.5(3) 

HBA 

6.0(3) 

13.7(3) 

P 

6.4(1) 

8.8(1) 

L 

6.7(3) 

8.7(3) 

HBA 

2(2) 

25(2) 

P 

6(1) 

19(1) 

L 

6.7(3) 

20.1(4) 

LBA 

14(2) 

20(2) 

HBA 

9.8(6) 

14.1(6) 

P 

7(2) 

11(1) 

L 

6.7(3) 

8.8(3) 

HBA 

4.6(3) 

22.9(3) 

P 

8.2(6) 

19.2(6) 

L 

4.2(3) 

10.9(3) 

LBA 

6.3(6) 

26.0(8) 

HBA 

7.0(6) 

20.4(6) 

P 

13(1) 

23(3) 

L 

7(2) 

13.4(4) 

LBA 

33.0(3) 

65(1) 

HBA 

16(5) 

50(5) 

P 

12.3(7) 

28.5(7) 

L 

11.6(3) 

29.2(4) 

LBA 

4.2(9) 

12(1) 

HBA 

4.2(6) 

9.8(6) 

P 

2.9(6) 

7.1(6) 

L 

5.3(3) 

8.4(3) 

LBA 

14.1(3) 

18.3(3) 

HBA 

10.9(3) 

14.1(3) 

P 

9.8(1) 

12.7(1) 

L 

1.4(3) 

8.8(3) 

LBA 

1.4(6) 

22.5(6) 

HBA 

1.8(3) 

17.6(3) 

P 

12.4(7) 

14.6(7) 

L 

10.9(3) 

13.0(3) 

HBA 

1.4(6) 

9.8(6) 

L 

3.2(4) 

10(1) 

HBA 

4(1) 

8(1) 


2.71(3) 

2.54(9) 

-0.17(9) 

-0.12(6) 

5.5(3) 

3.98(9) 

3.62(9) 

-0.39(8) 

-0.47(6) 

9.7(7) 

3.5(3) 

2.6(2) 

1.66(9) 

-0.05(4) 

-0.16(2) 

4.3(9) 

3.81(9) 

2.44(3) 

2.41(9) 

0.0(3) 

-0.06(8) 

7.0(7) 

5.5(3) 

5.6(1) 

-0.18(7) 

-0.21(5) 

5.7(7) 

3.9(1) 

3.3(4) 

2.44(9) 

-0.0(1) 

-0.31(3) 

6.35(9) 

5.3(1) 

3.03(9) 

0.3(1) 

-0.18(3) 

7.2(2) 

5.7(1) 

6.6(8) 

3.7(1) 

-0.31(2) 

-0.22(2) 

18.1(4) 

14(1) 

7.9(2) 

8.1(1) 

0.1(1) 

-0.08(6) 

3.5(4) 

2.7(1) 

2.0(1) 

2.35(9) 

-0.16(3) 

-0.17(2) 

5.08(9) 

3.91(9) 

3.52(3) 

2.44(9) 

0.2(1) 

-0.15(2) 

6.3(1) 

4.88(9) 

4.1(1) 

3.62(9) 

0.4(4) 

0.0(1) 

2.7(1) 

2.8(3) 

-(-) 

-(-) 

2.3(3) 

-0.0(1) 

-0.07(9) 


B1322+83 
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PSR Name 

Band 

W50 

[deg] 

i^io 

[deg] 

wwlP 

% 

^50 

<5io 


HBA 

11(1) 

22(1) 

6.3(3) 




P 

8(1) 

22(1) 

6.4(5) 




L 

10.2(5) 

19(1) 

5.3(3) 



B1508+55 





0.37(8) 

0.07(4) 


LBA 

scattered 






HBA 

4.2(3) 

9.5(3) 

2.64(9) 




P 

4.6(1) 

12.6(1) 

3.49(2) 




L 

8.1(3) 

12.3(3) 

3.42(9) 



B1530+27 





-0.42(9) 

-0.23(5) 


HBA 

12.0(6) 

20.4(8) 

5.7(2) 




P 

8(1) 

15.1(7) 

4.2(2) 




L 

4.6(3) 

12.0(3) 

3.3(1) 



B1540-06 





-0.2(1) 

-0.06(7) 


LBA 

7(1) 

18(3) 

5.1(8) 




HBA 

2(1) 

11(1) 

3.1(3) 




P 

3.3(2) 

8.3(2) 

2.31(5) 




L 

2.8(3) 

8.8(3) 

2.44(9) 



B1541+09 





0.52(2) 

0.44(4) 


LBA 

scattered 






HBA 

18.6(3) 

36.6(5) 

10.2(1) 




P 

17(3) 

71(3) 

19(1) 




L 

61(1) 

97(3) 

26(1) 



B1604-00 





-0.1(6) 

-0.04(4) 


LBA 

8(2) 

11(2) 

3.1(6) 




HBA 

9.5(3) 

14.1(3) 

3.91(9) 




P 

10.5(6) 

15.4(6) 

4.3(1) 




L 

7.4(3) 

12.7(3) 

3.52(9) 



B1612+07 





-0.0(4) 

-0.0(1) 


HBA 

1.4(6) 

7.0(6) 

2.0(1) 




P 

2.3(7) 

7.7(8) 

2.1(2) 




L 

1.8(3) 

6.7(4) 

1.9(1) 



J1627+1419 





-(-) 

-(-) 


HBA 

23(1) 

32(5) 

8(1) 



B1633+24 





-0.1(1) 

-0.5(7) 


LBA 

14(2) 

18(2) 

5.1(7) 




HBA 

8(2) 

28(2) 

7.8(7) 




P 

9(1) 

33(2) 

9.4(7) 




L 

9(1) 

20.1(4) 

5.6(1) 



B1642-03 





0.1(1) 

0.17(5) 


HBA 

2.8(3) 

7.7(3) 

2.15(9) 




P 

3.0(1) 

6.3(1) 

1.74(5) 




L 

3.5(3) 

10.2(3) 

2.83(9) 



J1645+1012 





-0.2(2) 

-0.1(1) 


HBA 

9(1) 

18(1) 

5.1(3) 




L 

6(1) 

16(2) 

4.5(6) 



J1652+2651 





-0.17(8) 

-0.16(7) 


HBA 

16(1) 

21(1) 

5.9(3) 




P 

14(1) 

16(1) 

4.7(2) 




L 

11.6(3) 

14.4(4) 

4.0(1) 



B1737+13 





0.2(2) 

-0.02(7) 


HBA 

5(1) 

19(1) 

5.5(3) 




P 

13.8(5) 

22(1) 

6.3(3) 




L 

13(1) 

19.7(7) 

5.5(2) 



B1749-28 





-0.3(1) 

-0.25(5) 


HBA 

6.0(3) 

12.0(3) 

3.32(9) 




P 

4.2(1) 

8.0(1) 

2.22(3) 




L 

3.5(3) 

7.4(3) 

2.05(9) 



B1818-04 





-(-) 

-(-) 


HBA 

scattered 






P 

5.9(7) 

16.3(7) 

4.5(1) 








A(&A-25196_ap, Online Material p 31 


Table B.2. Continued. 


PSR Name 

Band 

W50 

ww 

wwlP 

^50 

<5io 



[deg] 

[deg] 

% 




L 

6.0(3) 

11.3(3) 

3.13(9) 



B1821+05 

HBA 

3.5(3) 

8.8(3) 

2.44(9) 

0.95(9) 

0.49(4) 


P 

3.2(3) 

19.3(8) 

5.4(2) 




L 

22.5(3) 

28.9(5) 

8.0(1) 



B1822-09 (IP) 

LBA 

8(1) 

33(2) 

9.4(6) 

-0.2(1) 

0.4(1) 


HBA 

4(1) 

9(1) 

2.7(3) 




P 

4(1) 

9(1) 

2.6(5) 




L 

4.2(3) 

21.5(3) 

5.96(9) 



B1831-04 

HBA 

scattered 



-(-) 

-(-) 


P 

50(7) 

122(15) 

34(4) 




L 

108(2) 

132.0(9) 

36.7(2) 



B1839+09 

HBA 

8(2) 

19(2) 

5.5(7) 

0.3(2) 

-0.0(1) 


P 

3.9(7) 

13.2(8) 

3.7(2) 




L 

8.8(3) 

14.8(9) 

4.1(2) 



B1839+56 

LBA 

2.1(6) 

5.6(6) 

1.6(1) 

0.4(1) 

0.21(6) 


HBA 

1.8(3) 

3.5(3) 

0.98(9) 




P 

4.3(8) 

11.4(6) 

3.2(1) 




L 

6.3(3) 

9.9(3) 

2.74(9) 



B1842+14 

LBA 

scattered 



0.19(6) 

0.12(4) 


HBA 

33.8(3) 

61.9(3) 

17.19(9) 




P 

6.7(5) 

12.0(5) 

3.3(1) 




L 

7.7(3) 

12.4(3) 

3.44(9) 



B1848+12 

HBA 

4.2(6) 

11(1) 

3.3(3) 

-0.2(2) 

-0.3(1) 


P 

2.4(4) 

5.4(4) 

1.5(1) 




L 

2.8(3) 

6.0(5) 

1.7(1) 



B1848+13 

HBA 

8(2) 

19(2) 

5.5(7) 

-0.2(2) 

-0.2(1) 


P 

5(1) 

14(2) 

4.1(6) 




L 

4.9(3) 

12.0(4) 

3.3(1) 



B1857-26 

HBA 

scattered 



-(-) 

-(-) 


P 

19.3(9) 

39.6(9) 

11.0(2) 




L 

30.3(3) 

37.0(4) 

10.3(1) 



B1905+39 

HBA 

19.7(6) 

26.0(9) 

7.2(2) 

-0.11(4) 

-0.14(6) 


P 

16.5(7) 

20.1(7) 

5.6(1) 




L 

15.1(3) 

19(1) 

5.5(3) 



B1907+00 

HBA 

scattered 



-(-) 

-(-) 


P 

1.7(4) 

3.9(4) 

1.1(1) 




L 

1.8(3) 

8.8(7) 

2.4(1) 



B1907+02 

HBA 

scattered 



-(-) 

-(-) 


P 

2.7(8) 

9.9(8) 

2.7(2) 




L 

2.1(3) 

12(1) 

3.4(2) 



B1907+10 

HBA 

scattered 



-(-) 

-(-) 


P 

5.7(5) 

15.6(5) 

4.3(1) 




L 

6.0(3) 

16.9(3) 

4.69(9) 



B1911-04 

HBA 

scattered 



-(-) 

-(-) 


P 

2.95(9) 

5.71(9) 

1.59(2) 




L 

3.2(3) 

7.0(3) 

1.95(9) 
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Table B.2. Continued. 


PSR Name 

Band 

W50 

u'lo 

wwlP 

^50 

<5io 



[deg] 

[deg] 

% 



B1914+09 

HBA 

scattered 



-(-) 

-(-) 


P 

4(1) 

13(1) 

3.6(5) 




L 

12.0(3) 

16.9(9) 

4.7(2) 



B1915+13 

HBA 

scattered 



-(-) 

-(-) 


P 

4(1) 

14(1) 

4.1(3) 




L 

7.4(3) 

15.5(3) 

4.3(9) 



B1917+00 

HBA 

scattered 



-(-) 

-(-) 


P 

1.4(6) 

9.1(6) 

2.5(1) 




L 

1.4(3) 

10.2(3) 

2.83(9) 



B1918+26 

HBA 

1(1) 

9(1) 

2.7(3) 

-0.1(2) 

-0.2(1) 


P 

7.5(7) 

11(2) 

3.1(5) 




L 

6.0(6) 

7.4(3) 

2.1(1) 



B1919+21 

LBA 

5.3(3) 

11.3(3) 

3.13(9) 

0.11(5) 

-0.02(3) 


HBA 

7.0(3) 

10.2(3) 

2.83(9) 




P 

7.0(6) 

9.8(6) 

2.7(1) 




L 

8.4(3) 

10.6(3) 

2.93(9) 



B1920+21 

HBA 

scattered 



-(-) 

-(-) 


P 

2.7(6) 

13.0(6) 

3.6(1) 




L 

5.3(3) 

17.2(5) 

4.8(1) 



B1923+04 

HBA 

9.1(9) 

13(1) 

3.7(3) 

-0.3(1) 

-0.15(9) 


P 

4.9(7) 

10.7(7) 

3.0(2) 




L 

4.2(3) 

9.1(3) 

2.54(9) 



B1929+10(IP) 

LBA 

11(5) 

50(7) 

14(2) 

0.13(8) 

-0.18(3) 


HBA 

6.3(4) 

28.1(9) 

7.8(2) 




P 

9(1) 

21.8(3) 

6.06(9) 




L 

8.8(3) 

17.9(3) 

4.99(9) 



B1953+50 

LBA 

16(4) 

14(3) 

3(1) 

-0.2(1) 

-0.13(5) 


HBA 

3.2(3) 

10.9(3) 

3.03(9) 




P 

2(1) 

6(1) 

1.9(3) 




L 

3.2(3) 

8.1(3) 

2.25(9) 



B2021+51 

HBA 

5(1) 

25(1) 

7.0(3) 

0.3(1) 

-0.22(5) 


P 

4.5(5) 

21.2(5) 

5.9(1) 




L 

8.1(3) 

15.8(3) 

4.4(9) 



B2022+50 (IP) 

HBA 

5(2) 

22(3) 

6(1) 

-0.4(4) 

-0.5(1) 


P 

4(2) 

14(2) 

4.0(5) 




L 

2.5(3) 

8.1(3) 

2.25(9) 



B2043-04 

HBA 

4(1) 

9(1) 

2.5(4) 

-0.1(1) 

-0.0(1) 


P 

4.2(3) 

7.6(3) 

2.1(1) 




L 

3.9(3) 

7.7(3) 

2.15(9) 



B2044+15 

HBA 

2.1(6) 

9(1) 

2.7(3) 

0.1(3) 

-0.07(7) 


P 

2.8(9) 

16.9(6) 

4.7(1) 




L 

2.8(3) 

14.1(3) 

3.91(9) 



B2106+44 

HBA 

scattered 



-(-) 

-(-) 


P 

18(3) 

50(3) 

14.0(9) 




L 

21.1(3) 

31.7(3) 

8.8(9) 



B2110+27 





-0.3(2) 

-0.2(1) 
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PSR Name 

Band 

W50 

ww 

wwlP 

^50 

<5io 



[deg] 

[deg] 

% 




HBA 

3.5(6) 

6.3(6) 

1.8(1) 




P 

3.3(3) 

6.3(3) 

1.74(9) 




L 

1.8(3) 

4.6(3) 

1.27(9) 



B2113+14 

HBA 

scattered 



-(-) 

-(-) 


P 

6(1) 

14(1) 

3.9(5) 




L 

7.0(3) 

19(1) 

5.3(3) 



B2122+13 

HBA 

1(5) 

22(2) 

6.2(6) 

-0.2(2) 

-0.2(1) 


P 

15(2) 

16(2) 

4.7(6) 




L 

11.6(3) 

14(1) 

4.1(3) 



B2148+63 

HBA 

scattered 



-(-) 

-(-) 


P 

18(1) 

31(1) 

8.8(3) 




L 

11.3(3) 

19.0(3) 

5.28(9) 



B2154+40 

HBA 

7.4(3) 

33(1) 

9.4(3) 

-(-) 

-(-) 


P 

8.5(7) 

26.4(7) 

7.3(2) 




L 

14.8(3) 

23.9(3) 

6.65(9) 



B2217+47 

LBA 

scattered 



0.1(1) 

0.1(5) 


HBA 

4.6(3) 

9.1(3) 

2.54(9) 




P 

4.7(2) 

8.8(2) 

2.44(7) 




L 

5.3(3) 

11.3(3) 

3.13(9) 



B2224+65 

HBA 

11.6(3) 

22.9(9) 

6.3(2) 

0.48(3) 

0.22(4) 


P 

7(2) 

48(2) 

13.3(6) 




L 

34.5(3) 

44.0(4) 

12.2(1) 



B2227+61 

HBA 

scattered 



-(-) 

-(-) 


P 

20(1) 

26(1) 

7.4(3) 




L 

2.8(3) 

23(1) 

6.6(3) 



J2248-0101 

HBA 

scattered 



-(-) 

-(-) 


P 

7.4(8) 

13(1) 

3.6(3) 




L 

6.3(3) 

13(1) 

3.7(2) 



J2253+1516 

HBA 

2.1(7) 

22(2) 

6.3(6) 

-(-) 

-(-) 

B2255+58 

HBA 

scattered 



-(-) 

-(-) 


P 

6(1) 

24(1) 

6.7(4) 




L 

9.1(3) 

20.4(3) 

5.67(9) 



B2303+30 

LBA 

5(2) 

11(2) 

3.3(6) 

-0.2(1) 

-0.2(6) 


HBA 

5.6(3) 

12.3(5) 

3.4(1) 




P 

4.5(2) 

8.1(2) 

2.26(7) 




L 

3.2(3) 

7.4(3) 

2.05(9) 



B2306+55 

HBA 

4.9(3) 

30.6(4) 

8.5(1) 

0.43(7) 

-0.14(2) 


P 

20(1) 

24(1) 

6.8(4) 




L 

16.9(3) 

22.2(4) 

6.2(1) 



B2310+42 

HBA 

scattered 



0.07(5) 

-0.07(3) 


HBA 

8.4(3) 

18.3(3) 

5.08(9) 




P 

8(1) 

13(1) 

3.9(2) 




L 

9.9(3) 

15.5(3) 

4.3(9) 



B2315+21 

HBA 

5.6(6) 

7.7(6) 

2.1(1) 

-0.1(1) 

-0.0(1) 


P 

2(1) 

7(1) 

2.0(3) 




L 

4.2(3) 

7.0(3) 

1.96(9) 
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Table B.2. Continued. 


PSR Name 

Band 

W50 

[deg] 

i^io 

[deg] 

wwlP 

% 

dso 

<5io 

B2334+61 





0.5(1) 

0.0(1) 


HBA 

8(2) 

25(2) 

7.0(7) 




P 

8.1(8) 

22(2) 

6.2(7) 




L 

17.2(4) 

25(1) 

7.1(3) 
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Table B.3. Peak ratios. For pulsars with double and multiple peaks (Col. 1), for each frequency band (Col. 2) the amplitudes of the two most 
prominent peaks normalised to the amplitude of the main peak of the profile are indicate d (Cols. 3 and 4), ordered by spin phase (j) (Cols. 5 and 
6), and their r atio P2/P1 (Col. 7). The horizontal line in the table marks the change between previously studied cas es and new (as in Fig. [6l(. 


PSR Name 

Band 

PI 

P2 

01 

02 

P2/P1 

B0059+65 

HBA 

0.28 

0.96 

0.47 

0.52 

3.41 


P 

0.96 

0.64 

0.48 

0.52 

0.66 


L 

0.99 

0.75 

0.48 

0.52 

0.77 

B0301+19 

HBA 

1.00 

0.67 

0.48 

0.52 

0.68 


P 

0.89 

1.00 

0.48 

0.52 

1.12 


L 

0.66 

0.98 

0.49 

0.51 

1.49 

B0320+39 

LBA 

0.57 

0.68 

0.48 

0.50 

1.19 


HBA 

0.54 

0.79 

0.49 

0.50 

1.46 


L 

0.95 

0.32 

0.50 

0.51 

0.33 

B0329+54 

HBA 

1.00 

0.12 

0.50 

0.53 

0.12 


P 

1.00 

0.16 

0.50 

0.53 

0.17 


L 

0.99 

0.30 

0.50 

0.52 

0.30 

B0525+21 

LBA 

0.87 

0.70 

0.52 

0.58 

0.80 


HBA 

0.99 

0.75 

0.52 

0.58 

0.76 


P 

0.98 

0.92 

0.53 

0.57 

0.93 


L 

0.73 

1.01 

0.53 

0.57 

1.39 

B0834+06 

LBA 

0.97 

0.70 

0.50 

0.51 

0.72 


HBA 

1.01 

0.62 

0.50 

0.51 

0.62 


P 

1.00 

0.60 

0.50 

0.51 

0.60 


L 

0.93 

0.65 

0.49 

0.51 

0.70 

B0917+63 

LBA 

0.62 

0.85 

0.48 

0.52 

1.37 


HBA 

0.83 

0.95 

0.49 

0.51 

1.15 


P 

0.58 

1.00 

0.49 

0.51 

1.72 


L 

1.01 

0.93 

0.49 

0.50 

0.92 

B1133+16 

LBA 

0.71 

0.97 

0.49 

0.52 

1.35 


HBA 

0.98 

0.81 

0.49 

0.52 

0.82 


P 

0.98 

0.71 

0.49 

0.51 

0.72 


L 

0.99 

0.30 

0.49 

0.51 

0.31 

B1237+25 

LBA 

0.99 

0.34 

0.49 

0.53 

0.34 


HBA 

1.00 

0.44 

0.48 

0.52 

0.45 


P 

1.00 

0.61 

0.49 

0.52 

0.61 


L 

0.91 

0.98 

0.48 

0.51 

1.07 

B1604-00 

LBA 

1.10 

0.76 

0.49 

0.51 

0.69 


HBA 

0.52 

0.82 

0.49 

0.51 

1.57 


P 

0.73 

0.61 

0.49 

0.51 

0.83 


L 

0.56 

0.79 

0.49 

0.51 

1.42 

B2044+15 

HBA 

0.45 

0.97 

0.47 

0.51 

2.18 


P 

0.37 

1.00 

0.48 

0.51 

2.70 


L 

0.34 

0.99 

0.49 

0.51 

2.95 

B0402+61 

HBA 

0.66 

0.71 

0.49 

0.52 

1.08 


P 

0.86 

0.44 

0.49 

0.51 

0.51 


L 

0.48 

0.71 

0.51 

0.52 

1.50 

B0447-12 

HBA 

0.95 

0.43 

0.48 

0.50 

0.45 


P 

0.81 

0.75 

0.49 

0.52 

0.92 


L 

0.70 

0.87 

0.48 

0.51 

1.24 

B0450+55 

LBA 

1.02 

0.22 

0.47 

0.53 

0.22 


HBA 

0.89 

0.67 

0.49 

0.55 

0.75 


P 

0.99 

0.24 

0.49 

0.53 

0.24 


L 

0.48 

0.60 

0.46 

0.50 

1.27 

B0523+11 

HBA 

0.94 

0.84 

0.48 

0.52 

0.89 


P 

1.00 

0.67 

0.48 

0.52 

0.67 


L 

0.61 

0.98 

0.48 

0.51 

1.60 

B0609+37 

HBA 

0.36 

1.01 

0.48 

0.50 

2.77 


P 

0.26 

0.76 

0.49 

0.51 

2.91 


L 

0.32 

0.63 

0.50 

0.51 

2.00 

B0626+24 

HBA 

0.89 

0.27 

0.48 

0.50 

0.31 


P 

0.95 

0.41 

0.49 

0.50 

0.43 
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Table B.3. Continued. 


PSR Name 

Band 

PI 

P2 

01 

02 

P2/P1 


L 

0.48 

0.75 

0.49 

0.50 

1.57 

B0643+80 

HBA 

0.83 

0.95 

0.49 

0.50 

1.14 


P 

0.95 

0.48 

0.49 

0.50 

0.50 

B0906-17 

HBA 

0.95 

0.19 

0.49 

0.51 

0.20 


P 

0.95 

0.21 

0.49 

0.51 

0.22 


L 

0.32 

0.95 

0.48 

0.49 

3.00 

B0943+10 

LBA 

0.95 

0.24 

0.50 

0.54 

0.25 


HBA 

0.95 

0.41 

0.50 

0.52 

0.43 


P 

0.91 

0.76 

0.49 

0.51 

0.83 

B1530+27 

HBA 

0.83 

0.74 

0.49 

0.52 

0.90 


P 

1.02 

0.48 

0.50 

0.51 

0.47 

J1652+2651 

HBA 

0.89 

1.00 

0.48 

0.52 

1.13 


P 

0.82 

1.01 

0.48 

0.51 

1.23 


L 

0.65 

0.84 

0.49 

0.51 

1.30 

B1737+13 

HBA 

0.21 

1.03 

0.47 

0.51 

4.88 


P 

0.63 

1.01 

0.48 

0.51 

1.61 


L 

0.93 

0.57 

0.48 

0.51 

0.62 

B1905+39 

HBA 

0.71 

0.38 

0.48 

0.49 

0.53 


P 

0.95 

0.60 

0.48 

0.51 

0.63 


L 

0.95 

1.00 

0.49 

0.51 

1.05 

B1918+26 

HBA 

1.00 

0.26 

0.50 

0.51 

0.26 


P 

1.07 

1.03 

0.49 

0.51 

0.96 


L 

0.92 

0.72 

0.50 

0.51 

0.78 

B1919+21 

LBA 

0.57 

0.59 

0.49 

0.50 

1.03 


HBA 

0.95 

0.71 

0.49 

0.51 

0.75 


P 

0.95 

0.59 

0.49 

0.50 

0.62 


L 

0.74 

0.95 

0.50 

0.51 

1.29 

B2122+13 

HBA 

0.40 

0.99 

0.47 

0.51 

2.47 


P 

0.89 

0.96 

0.48 

0.51 

1.09 


L 

0.74 

0.91 

0.48 

0.51 

1.23 

B2154+40 

HBA 

0.78 

0.35 

0.50 

0.53 

0.45 


P 

0.89 

0.46 

0.49 

0.52 

0.52 


L 

0.62 

0.69 

0.49 

0.50 

1.12 

B2306+55 

HBA 

0.46 

1.00 

0.46 

0.52 

2.18 


P 

0.79 

0.98 

0.47 

0.52 

1.24 


L 

0.97 

0.87 

0.49 

0.52 

0.89 

B2315+21 

HBA 

0.85 

0.67 

0.50 

0.51 

0.79 


P 

0.32 

1.11 

0.49 

0.50 

3.50 


L 

0.38 

0.95 

0.49 

0.50 

2.50 






